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Offner spectrometer of VIRTIS-M: 1) mirror; 2) diffrac¬ 
tion grating; 3) the entrance slit and electromechan¬ 
ical shutter; 4) infrared focal plane; 5) visible focal 
plane; 6) position of the mirror on visible focal plane 
(not shown); 7) mechanical structure (Selex ES). . . . 
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Equivalent to hgure 12.31 for one single sample. The 
rise above 4 pm is due to the thermal background. 
The different electronic offset along the bands caused 
by the multiplexer readout circuit can also be noted. 
Left panel; Mean Dark current versus integration 
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first 250 bands. The solid line is a linear best fit. The 
Pearson correlation coefficient is used for the analy¬ 
sis. The intercept of the line represents the mean 
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points of left panel as a function of the Focal Plane 

Temperature. 

This dark frame of VIS channel clearly shows an ex¬ 
ternal source of light at longer wavelengths. Frames 
like this are not taken into account in the analysis. 
Left panel : Mean Dark current versus integration 
time. The solid line is a linear best ht. The Pearson 
correlation coefficient is used for the analysis. The 
intercept of the line represents the mean electronic 
offset in DN while the slope corresponds to the dark 
current rate in DN/s. Right panel: same points of left 
panel as a function of the Focal Plane Temperature. 
The negative correlation observed is due to the low 
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2.8 Electronic offset for even bands (IR channel) in func¬ 

tion of the temperature of the focal plane. There is 
not signihcant correlation. 

2.9 Electronic offset for odd bands (IR channel) in func¬ 

tion of the temperature of the focal plane. There is a 
negative correlation. 

2.10 Mean offset of VIS channel vs Temperature of the 

focal plane. The correlation is negative. 

2.11 The simulated DN as a result of four different temper¬ 
atures of the spectrometer, with an integration time 
of 1 s. Thermal Background affects the spectrum 
at decreasing wavelength for increasing temperature. 
The pixels are saturated where the line is missing. 
This simulation is performed for the central sample 
of the slit. The result is similar for the other samples 

2.12 The simulated S/N (blue line) is over plotted to the 

S/N extracted from the visible lamp signal (black 
line) acquired during internal calibration mode with 
IR channel. The integration time is 0.1 s. The av¬ 
erage overestimation of the simulator is 30%. This 
comparison also represents a satisfactory test for the 
saturation of IR channel (where the S/N is set to 0). 
The band at ~ 2800 nm is the limit before which the 
signal is saturated. This is what emerges from the 
real data and the simulation. 

2.13 The simulated S/N (blue line) is over plotted to the 

S/N extracted from the infrared lamp signal (black 
line) acquired during internal calibration mode with 
IR channel. The average overestimation of the simu¬ 
lator is 40%. Integration time = 0.5 s. 

2.14 The simulated S/N (blue line) is overplotted to the 
S/N extracted from the Lutetia spectral cube (black 
line) in IR channel. The integration time is 0.7 s. 
The average overestimation of the simulator is 40%. , 

2.15 The simulated S/N (blue line) is over plotted to the 

S/N extracted from the visible lamp signal (black 
line) acquired during internal calibration mode with 
VIS channel. The integration time is 1 s. The average 
overestimation of the simulator is 15%. 
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2.16 The simulated S/N (blue line) is over plotted to the 
S/N extracted from the infrared lamp signal (black 
line) acquired during internal calibration mode with 
VIS channel. The integration time is 20 s. The aver¬ 
age overestimation of the simulator is 10%. 


2.17 The simulated S/N (blue line) is over plotted to the 
S/N extracted from the Lutetia spectral cube (black 
line) in VIS channel. The integration time is 1 s. 
The simulator underestimates the S/N in average by 
a factor 15%. This underestimation of S/N is due to 
the underestimation of the signal. In fact the VIS 
channel is affected by c ontamination frorn the high 
orders of the IR channel f Filacchiond . 20061) . and this 
contamination is not well characterized by the ITF. 
This is a further source of error that is not taken into 
account in this analysis. 
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The dehned SSA of the Dark Terrain is plotted in 
black. Two different mixing are over plotted. The 
Dark Terrain is mixed with water ice, carbon dioxide 
and methanol, with different abundances. 

Retrieved radiance (up panel) and reflectance (down 
panel) with error bars. In this example the simu¬ 
lated composition is: intimate mixing with 5% water 
ice, 5% carbon dioxide, 5% methanol, 85% dark ter¬ 
rain. The heliocentric distance is hxed to 2.5 AU, 
from which the resulting temperature of the surface 
is 189 K. The integration time is fixed to 2 s. The 
error bars are calculated as explained in section 4.3. 
The error bars for radiance spectrum are enhanced 
by a factor 10 to make them clearly visible. Other 
parameters are fixed as stated in the text. 
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Relative error on the retrieved band area versus dif¬ 
ferent integration times, for heliocentric distance of 
3.2 AU. The four panels represent different mixing 
modalities: intimate (up), areal (down), 1% abun¬ 
dances (left), 5% abundances (right). 

Same of Figure [231 for heliocentric distance of 2.5 AU. 
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3.5 Same of Figure 13.31 for heliocentric distance of 1.5 
AU. The missing points represents saturated pixels in 
correspondence of the absorption bands of the end- 
members. 
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4.1 


As a first guess the reflectance spectrum is modeled 
taking into account the dark terrain and the optical 
constants of the water ice and the carbon dioxide. 
Up: the synthetic spectrum and the resulting best 
fit (clear line). The reduced chi-square results to be 
3.6. Down: the residuals spectrum smoothed with 
a 10 bands running box. Only the residuals having 
the signal overestimated by the model are shown to 
enlighten absorption bands we are missing. We con¬ 
sider significant a residual if it overcome 2 a. We can 
notice a significant residual at ~ 3.5 /xm. We can 
argue we are missing the absorption band of methanol 
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4.2 We repeat the procedure of Figure ITT] including methanol. 
The resulting reduced chi-square is 1.1. The bottom 
panel show the residuals spectrum keeping the same 
range of Figure 14.11 for the sake of clarity. There 
are no significant residuals. We can argue we have 
achieved a good model. 
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4.3 Four thousand noisy spectra are simulated, and the 
plots show the probability density function (80 bins) 
of the retrieved abundance of water ice, being the 
value of the model in input equal to 5%. The black 
line shows the resulting distribution in the case is 
defined according to Eq. 14.11 The clear line shows the 
resulting distribution in the case the a in Eq. 14.11 is 
not coming from the S/N simulator but the observed 
data is used as weight. The theoretical gaussian dis¬ 
tributions over plotted (dashed line) have the same 
mean and dispersion of the resulting distributions. 
The position of the bars are in the points of the the¬ 
oretical distributions such that they contain the 95% 
of the probability distribution. The distance of the 
bars from the mean value is indicated in the figure. 
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4.4 Schematic diagram of tilt coordinates on an element 
of surface S. Being N the normal to the surface, Nt 
the normal tilted, i and e the incidence and emission 
angle, it and ct the incidence and emission angle 
after the tilt of the normal, ip and p' respectively the 
azimuth of the tilt and the azimuth of the emission 
direction from the direction of the incidence angle; 9 
is the tilt angle, and g is the phase angle. 
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Deep Impact spacecraft. 

Comparison of a small part of (left) Tempel 1 with 
(right) Hartley 2 at approximately the same image 
scale and with nearly identical instruments. Left: Im- 
pactor Targeting Sensor (ITS) image iv9000675, 9.1 
m/pixel. Right: M RI image mv5004032, 8.5 m/pixel. 

Sun is to the right f A’Hearn et ahl. I201D . 

The measured radiance (black line) is modeled in 
terms of reflected sunlight and thermal emission. . . . 
Left panel: HRIV acquisition. Exposure ID = 9000901, 
right panel: same image after lowering the spatial res¬ 
olution to match that of the HRII channel. 

Average spectra of nucleus regions shown in Fig. 15.61 
Panel A: image of the nucleus of the comet at A = 
0.35 /rm (exposure ID 9000906), rescaled to the same 
resolution of the hyperspectral cubes. Panel B: spec¬ 
tral classihcation of the comet’s nucleus: in blue the 
icy regions, in green the low reddening regions, in or¬ 
ange the high reddening regions, in grey the regions 
which presents unfavourable angles for the analysis. . 
Intimate (left) and areal (right) mixing results: icy 
region (corresponding to pixels marked in blue in Fig. 
15.61) average observed spectrum (black curve); areal 
mixing spectra with 2.4, 1.2 and 0.4% of water ice are 
indicated by red, green and blue curves, respectively. 
The best-ht are obtained respectively with 1.2% and 

0.4% water ice-rich spectrum. 

Water ice abundance map, in case we model the mix¬ 
ing between Dark Terrain and pure water ice to be 
intimate. 
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5.9 Upper left panel; temperature T{K)] upper right 
panel: effective emissivity eg//; lower right panel: self 
heating; lower left panel: reddening (%/kA), it is cal¬ 
culated with respect the value of the reflectance at 
A = 1.7 um for a direct comparison with the work of 


Davidsson et al.l (120091) . l91 


5.10 Mean (left) and standard deviation (right) map of the 
temperature on the surface, according to the Gaus¬ 
sian distribution of Planck functions model. 93 


5.11 Water ice abundance and grain size map of Hartley 
2, in case we model the mixing between amorphous 
carbon and pure water ice to be intimate. 
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5.12 Upper panel: the measured spectrum and the re¬ 
trieved model. Lower panel: residuals analysis. The 
spectrum of residual is smoothed by a running box. 
A signihcant residual (> 2a) at A ^2.7 um is re¬ 


lated with the water vapour emission flA’Hearn et ah 


201 ll) which is not taken into account in the model 


The best £t is performed skipping the spectral range 
2.55 - 2.8 fim (dashed spectrum), related to this wa¬ 
ter vapour emission. Retrieved water ice abundance: 
1 %; retrieved grain sizes: water ice = 100pm, amor¬ 
phous carbon = 4 pm. 
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5.13 Upper panel: the measured spectrum and the re¬ 
trieved model. Lower panel: residuals analysis. The 
spectrum of residual is smoothed by a running box. 
A signihcant residual (> 2a) at A ^2.7 um is re¬ 


lated with the water vapour emission flA’Hearn et ah 


201 ll) which is not taken into account in the model 


The best £t is performed skipping the spectral range 
2.55 - 2.8 pm (dashed spectrum), related to this water 
vapour emission. Retrieved water ice abundance: 8 
%; retrieved grain sizes: water ice = 7pm, amorphous 
carbon = 2 pm. 






















12 


LIST OF FIGURES 


5.14 The single scattering albedo (a;) calculated for the 
two comets. The high level of uncertainty coming 
from the thermal emission removal prevents the cal¬ 
culation of the SSA beyond 3.5 /rm. The SSA derived 
from Hartley 2 data presents contamination of water 
vapour emission (A ~ 2.7) not removed here. Super¬ 
imposed is the linear interpolation of the Tempel 1 
SSA, which has been used for the spectra simulation 
in Chapter [3l The two slopes are related to the two 
different channels of the instrument: HRIV and HRII 
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5.15 Upper left panel: temperature T{K); upper right 
panel: effective emissivity Ce//; lower right panel: self 
heating; lower left panel: reddening (%/kA), it is cal¬ 
culated with respect the value of the reflectance at 
A = 1.7 um for a direct comparison with the work of 


Davidsson et al.l (120091) 


5.16 Mean (left) and standard deviation (right) map of the 
temperature on the surface, according to the Gaus¬ 
sian distribution of Planck functions model (see sec¬ 
tion [51131 ) . 
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6.1 


Lutetia mean reflectance spectrum in the 1.0 - 5.0 
fim spectral range. Where the line is thin there are 
instrumental artifacts in the spectrum; the rise above 
4 fim is due to the surface thermal emission. In light 
color are shown the la spectra. The panel inside 
the plot shows a detail of the reflectance within the 
wavelength range taken into account for the analysis. 
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6.2 


The four hyperspectral cubes analyzed at A = 2.0/im. 

They are (from up to down) the four listed in table I6.1lll04 
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6.3 


6.4 


6.5 

6.6 
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Left panel: reflectance centered in 1.142 /im, filter 
0.1 fim wide, plotted against phase angle for all se¬ 
lected spectra. Scatter of the points is dne to different 
observation geometry (incidence and emission angle). 

In the zoom is shown a detail for small phase angles, 
three of the fonr hyperspectral cnbes cover this range. 

From 4° to 7.5° the points are the result of the aver¬ 
age over the whole disk (14 scans). Right panel: the 
best £t of the model. It is able to account for dif¬ 
ferent observation geometry. Most of the discrepancy 
between data and model are probably due to errors 

in the geometry information. Il06 

The retrieved parameters are plotted as a function 
of the wavelength. The mean slope 9 is expressed in 
radians, like the hs parameter. The b parameter is 
plotted with an offset of -0.1 for the sake of clarity. 

Error bars are calculated as discussed in section 14.21 
The map of artifact in the IR channel. It is the me¬ 
dian of the Eq. l6.2l for a large numbers of acquisitions. 
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The crack of the hlter is clearly visible. IllO 

The map of artifact in the VIS channel. It is the me¬ 


dian of the Eq. 16.21 for a large numbers of acquisitions. lllO 
An example of cleaned spectrum of Lutetia after de- 


spiking and artifact removal. Illl 

Upper panel: the measured spectrum and the re¬ 
trieved model. Lower panel: residuals analysis. The 
spectrum of residual is smoothed by a running box. 
Significant residuals (> 2a) affect many parts of the 
spectrum, which is in fact poorly htted by the model. 

For a complete discussion refer to the text. Retrieved 
abundances: tholin = 7.5 %, amorphous carbon = 

92.5 %; retrieved grain sizes: tholin = 2.5 pm, amor¬ 
phous carbon = 1.1 pm. Ill3 


7.1 Rosetta navigation camera image taken on 19 August 
2014 at about 79 km from comet 67P/Churyumov- 
Gerasimenko. The comet nucleus is about 4 km across. 

{http://www.esa.int/spaceinimages/Images/ 2 Oi 4 /O 8 /ConwEon_l 9-AugusE2014- 
_NavCam) . 118 
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Introduction 


Space science has certainly brought important breakthroughs in re¬ 
cent science history. It allowed to overcome the shield of terrestrial 
atmosphere to observe photons emitted by of planets, stars, galaxy 
and beyond, across the entire electromagnetic spectrum. In par¬ 
ticular for Solar System studies, it represents something very novel 
for astrophysics: the possibility to enter directly in contact, or very 
close, to the object under investigation, thanks to missions which 
are intergenerational in most cases, but whose travel time certainly 
worth the potential of discoveries. 

Not much time has passed since the hrst scientihc space mission 
has sent to Earth the hrst data. We can say that space exploration is 
still at the beginning, and every new data from in situ observations 
can bring a new comprehension of the Solar System, and thus of the 
part of the universe closer to us. 

Topic of this work are comets, small and elusive objects that may 
hold great secrets about the origin of the Solar System and life on 
Earth, being among the most primitive objects. 

The method of investigation addressed in this work is the visible- 
infrared spectrophotometry by remote sensing instruments. 

Spectrophotometry has already had very successful results in 
planetary exploration, in particular thanks to a new class of in¬ 
struments: imaging spectrometers, designed for the observation of 
remote planetary atmospheres and surfaces, capable to acquire hy- 
perspectral data with high spatial and spectral resolution. 

This work has to be intended as a preparation for the huge 
amount of data coming from VIRTIS (Visible and InfraRed Ther¬ 
mal Imaging Spectrometer) onboard Rosetta spacecraft, which will 
orbit and follows comet 67P/Churyumov-Gerasimenko (CG) during 
its journey toward the Sun, and will deliver a lander on the surface 
for the hrst time in space exploration history. 
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The context under which this mission moves its steps is described 
in the hrst chapter. 

In the second chapter the performances of the VIRTS instrument 
are analyzed in detail. In particular the modeling of the signal to 
noise ratio is the main argument of this chapter. 

The third chapter try to answer the question: how does VIRTIS 
manage cometary spectra? Here simulations of possible spectra of 
the comet’s nucleus are performed, which are useful for both a com¬ 
parison with real spectra, and for a planning of the observations. 
In this chapter it is also introduced the Hapke’s radiative transfer 
model we intend to use to invert acquired data to infer physical 
properties. 

The fourth chapter introduces a method for spectral modeling. 
It takes advantage of fast algorithms for optimization of least square 
problems to perform accurate retrieval of the investigated parame¬ 
ters. Moreover, the method includes the information on the instru¬ 
mental noise, permitting the analysis of the goodness of the models, 
and an estimation of the error of the retrieved parameters. 

In order to be prepared to manage VIRTIS data of the comet, we 
followed two different approaches: analysis of comet’s data acquired 
by another spectrometer and of VIRTIS data of another target (hfth 
and sixth chapters). 

In the hfth chapter data coming from previous missions to comets 
are analyzed: Deep Impact to Tempel 1 and its extended investi¬ 
gation to Hartley 2. Among the scientihc payload of Deep Imapct, 
the probe carries a spectrometer similar to VIRTIS (HRI). Its data 
allowed us to test the reliability of the method developed and to 
produce further investigations on those targets. 

In the sixth chapter analysis of another small body is presented: 
Lutetia asteroid. Its data are coming from VIRTIS itself thanks to 
a hyby of Rosetta during its cruise phase. 

This work have paved the way to the analysis of the hnal tar¬ 
get. The tools presented are currently used by the VIRTIS-Team to 
produce works on the comet, that are recommended to the reader. 
Since a complete analysis on the comet is outside the scope of this 
work, just preliminary results are shown here. 



Chapter 1 

Comets and space missions 


1.1 Cometary science history 


Comets are impressive celestial phenomena: bright blotches of light 
with long, beantifnl tails, snddenly appearing across the sky. They 
travel on elliptical orbits that swing aronnd onr Snn and retnrn 
back, in most cases beyond the orbit of Neptune. These incredibly 
elongated orbits mean that it can take years, from a few to hundreds 
or even thousands, for a comet to reappear to an observer on the 
Earth. 

The seemingly transient and unpredictable behavior of comets 
rendered them completely mysterious to ancient observer. Rather 
than being viewed as beautiful and intriguing, comets were instead 
regarded as omens of impending doom, or meteorological phenom¬ 
ena in terrestrial atmosphere and were not yet clearly established as 
astronomical bodies. 

Humans have been captivated by comets for much of history. 
Some archaeologists suggest that prehistoric rock paintings, found 
in several sites across the globe, may portray comets. Early rock 
carvings resembling comets have been found in Scotland dating back 
to the second millennium BC, and a comet-like shape found on rock 
carvings in Val Camonica, Italy, dates back to the late Iron Age (see 
Fig. dH]). 

The word “comet” comes from the Greek “kometes”, which liter¬ 
ally means “long-haired”, but the earliest extant reco rds o f corn et ary 
observations date from around 1000 BC in China (jHol. 19621) and 
probably from about the same time in Chaldea. 
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CHAPTER 1. COMETS AND SPACE MISSIONS 



Figure 1.1: Comet symbol on Rock 35 from Area di Foppe, Nadro di Ceto (Bres¬ 
cia, Italy). Riserva naturale Incision! rupestri di Ceto, Cimbergo e Paspardo. 
Copyright: Centro Camuno di Studi Preistorici, Dipartimento Valcamonica e 
Lombardia 

Ancient Greek philosophers put forward a number of different 
ideas to explain the physical nature of comets. To some, comets 
were celestial bodies just like planets. To others, they were linked 
to celestial mechanics such as planetary conjunctions, and to yet 
others they were burning clouds or optical phenomena in the Earth’s 
atmosphere. Although far from correct, the latter view was adopted 
by Aristotle in the fourth century BC. 

Aristotle suggested that they were “windy exhalations” from the 
Earth that reached out into our atmosphere. The Aristotelian ideas 
about planets and comets were upheld for an entire millennium dur¬ 
ing which there was little scientihc advancement in the held of as¬ 
tronomy. 

The only voice that dared contradict this view was that of Seneca, 
in the hrst century BC. According to Seneca’s Naturales Quaes- 
tiones, comets were more like planets. He also recognized that their 
sporadic behavior would make studying them quite complex, but 
hoped that future observations and investigations would shed light 
on their nature. 

Aristotle’s view of the heavens dominated almost undisputed 
throughout the Middle Ages. This situation started to change, 
however, with the dawn of the Renaissance in Europe. New ideas 
emerged as astronomers assembled an increasingly larger body of 
data. 




1.1. COMETARY SCIENCE HISTORY 
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In the fifteenth century, Italian scientist Paolo dal Pozzo Toscanelli 
conducted the first systematic observations of several comets. He 
recognized the importance of measuring both the orientation of the 
comet’s coma, the nebulous envelope surrounding the comet’s nu¬ 
cleus, and the position of the nucleus. 

In the sixteenth century, German astronomer Peter Apian and 
the Italian scientist Girolamo Fracastoro realized that a comet’s tail 
always points away from the Sun, a revelation that would be key to 
correctly interpreting these objects and their physics. 

The exceptionally accurate observations of the Danish astronomer 
Tycho Brahe (1578) initiated a new era for observational astronomy, 
as he demonstrated that the horizontal parallax of the bright Gomet 
G/1577 VI was certainly smaller than 15 arcmin, corresponding to 
a distance in excess of 230 Earth radii, consequently farther away 
than the Moon. 

The German astronomer Johannes Kepler observed several comets 
throughout his life, and tried to study their motion. However, he be¬ 
lieved comets to be interstellar objects moving along straight lines, 
and so never fully grasped their dynamics. 

Some progress on understanding cometary orbits came from Pol¬ 
ish astronomer Johannes Hevelius, who suggested in his Cometographia 
treatise, published in 1668, that comets moved on open-ended parabolic 
orbits, meaning that they would be seen once and never again. 

The English physicist and mathematician Isaac Newton, in his 
Philosophiae Naturalis Principia Mathematica, developed the bril¬ 
liant tool that could link all these observations, together with the 
studies of Galileo Galilei. With his new theory, Newton calculated 
the orbit of the Great Gomet of 1680 (Fig. II.2p . His results were 
in excellent agreement with observations. Shortly after, the English 
astronomer Edmond Halley applied Newton’s theory to the study 
of comets. This work would produce two major results: a strik¬ 
ing confirmation of Newton’s theory of universal gravitation, and 
a revolutionary change in our understanding of comets. With this, 
Halley was the first to suggest that comets may be periodic, moving 
along very elongated ellipses rather than parabolic paths. With an 
estimated period of about 76 years, Halley expected the comet to 
be visible again between the years of 1758 and 1759. Halley died 
before he could view it with his own eyes, but others were able to 
test and confirm his conjecture and, with it, Newton’s theory of 
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Figure 1.2: The Great Comet of 1680 over the Dutch town of Alkmaar, in a 
painting by Lambert Doomer. Rijksmuseum, Amsterdam. In public domain. 


gravity, demonstrating that it can be applied not only to planets 
but to other celestial bodies too. It also marked one of the greatest 
triumphs of science. 

Although it seemed like our understanding of comets was becom¬ 
ing more complete, one vital piece was missing from this puzzle: 
where do they come from? The debate shifted from nature to origin 
thanks to the German philosopher Immanuel Kant, who wrote his 
General History of Nature and Theory of the Heavens in 1755. In 
this early work, Kant suggested that the Sun and its planets formed 
from an extended diffuse nebula. However, the idea of an interstellar 
origin of comets prevailed until the second half of the nineteenth cen¬ 
tury, when the Italian astronomer Giovanni Schiaparelli suggested 
that comets belong to the Solar System and surround the Sun in an 
almost uniform cloud. He also established the link between comets 
and meteors, observing that the Perseid and Leonid meteor streams 
coincided with the orbits of Gomets 109P/Swift-Tuttle (1862 HI) 
and 55P/Tempel-Tuttle (1866 I) respectively. This was the proof 
that comets were indeed losing solid particles. 

After Halley, the German astronomer Johann Franz Encke (1820) 
was the second to successfully predict the return of a comet (in 
1822). Gomel 2P/ Encke has the shortest period of all known 
comets, 3.3 years, which provides similar Sun-Earth-comet conhgu- 
rations every 10 years. The comet was subsequently found to arrive 
systematically at perihelion about 0.1 days earlier than predicted, 
even when taking planetary perturbations into account. 




1.1. COMETARY SCIENCE HISTORY 


23 


Inspired by his observations of an asymmetric distribution of lu¬ 
minous matter in the head of Comet Halley in 1835, the German as¬ 
tronomer Friedrich Bessel (1836) interpreted this as a Sun-oriented 
asymmetric outflow and suggested that a non-gravitational effect 
might arise due to the rocket-type impulse imparted by such an 
outflow, possibly explaining perihelion shifts such as those observed 
for Comet 2P/Encke. It would take over a century for this idea to 
be fully accepted by the scientihc community because the existence 
of a solid body at the center of the coma was far from being unan¬ 
imously accepted. In fact, the theory that comets were a swarm of 
solid particles was the most favored by scientists at that time. 

While astronomers were making great progress using telescopes, 
the arrival of photography in the mid-nineteenth century opened up 
a new way to study our skies. The hrst comet to be photographed 
was the Great Comet of 1858, also known as Donati’ s Comet af¬ 
ter Italian astronomer Giovanni Donati who discovered it. A few 
years later, Donati was also the hrst astronomer to use spectroscopy 
to study the composition of a comet. By splitting the light from 
celestial bodies into its constituent colours through a prism, spec¬ 
troscopy allows astronomers to investigate the chemical composition 
of distant and otherwise inaccessible objects. Donati recorded the 
spectrum of a comet, now known as Comet C/1864 Nl, that had 
been discovered by German astronomer Ernst Wilhelm Tempel in 
1864. The spectrum contained three features that are now known 
to be produced by molecules of diatomic carbon {C 2 ) in the comet’s 
coma. Further observations around the turn of the twentieth cen¬ 
tury uncovered more about the chemical makeup of cometary co¬ 
mas, identifying sodium ions and a variety of carbon-, oxygen- and 
nitrogen-based molecules. 

In the hrst half of the twentieth century astronomers were col¬ 
lecting and studying more high-quality astronomical data than ever 
before, building up an impressive database. This allowed them to 
delve into the physical nature and origin of comets in great detail. 

In 1950, American astronomer Fred Whipple proposed a new 
model to describe comets. Rather than a loose collection of dust 
and debris kept together by ice, he suggested that comets have an 
icy nucleus, consisting primarily of frozen volatiles like water, carbon 
dioxide, methane, and ammonia, and containing only traces of dust 
and rock. 
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Then, the German astronomer Ludwig Biermann (1951) gave the 
correct explanation for the motions of features in cometary plasma 
tails caused by their interactions with a flow of charged particles 
emanating from the Sun’s surface (i.e., the solar wind). 

Finally, from dynamical studies of the distribution of semimajor 
axes of comets, came the identihcation by the Dutch astronomer 
Jan Oort (1950) of a distant population of comets now known as 
the Oort cloud, a huge cloud of “dormant” comets extending over 
a thousand times farther than the orbits of Neptune and Pluto. 
This cloud is not the only reservoir of comets in the Solar System; 
as suggested by Dutch astronomer Gerard Kuiper in 1951, most 
comets with a relatively short period are located in a flattened, 
ring-like distribution that begins just outside Neptune’s orbit, now 
known as the Kuiper belt. 

None of these ideas resulted directly from new observational ev¬ 
idence, but this was the hrst time that the known facts were ef¬ 
fectively combined, leading to a comprehensive description. A new 
picture of comets, and the existence of a family of celestial bodies, 
were suddenly revealed at the same time. 

Advances in technology towards the end of the twentieth century 
provided exciting new opportunities for comet scientists. Whereas 
before scientists had relied solely on ground-based telescopes to ob¬ 
serve and study comets, with the advent of the space age there was 
the novel opportunity to approach these icy bodies as they journeyed 
towards the inner regions of the Solar System, passing relatively 
close to Earth. 

In 1985, NASA’s International Gometary Explorer (IGE) became 
the hrst space mission to approach a comet passing approximately 
7800 kilometres from the nucleus of Gomet 21P/Giacobini-Zinner. 

Just one year later, when Halley’s Gomet returned to the vicinity 
of the Sun, an armada of spacecraft was sent to study the comet up 
close: these included two probes from Russia (Vega-1 and Vega-2), 
two from Japan (Sakigake and Suisei), and ESA’s Giotto spacecraft. 

Named after the Italian painter who, in 1303, depicted the star 
of Bethlehem as a comet, the Giotto mission obtained the closest 
pictures ever taken of a comet (Fig. II.3p . After Vega-1 and Vega-2 
came within 8900 and 8000 kilometres respectively from the nucleus 
of Halley’s Gomet in early March 1986, ESA scientists and engi¬ 
neers used the telemetry data from these probes to guide the Giotto 
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Figure 1.3: The nucleus of Comet Halley by ESA’s Giotto in 1986, from a 
distance of about 2000 km. Credit: ESA/MPS 


spacecraft even closer (reported from http://sci.esa.int/rosetta/). 

On 13 March 1986, Giotto flew past the comet’s nucleus at a 
distance of less than 600 kilometres, relaying unprecedented im¬ 
ages. It is only with this hrst spacecraft flyby that the idea of a 
single solid body at the center of the comet activity was widely ac¬ 


cepted by scientihc co mmunity 1 Keller et ah . 19861: Reinhard . 1986 


Sagdeev et al.l . ll986aHbl ) 


The nucleus was found to be larger (equivalent radius 5.5 km) 
and darker (albedo ~ 4% ) than expected. In the Giotto images, 
surface features (craters, ridges, ni ountains, etc.) and source regions 
were observed ( Keller et al.l. 1988h . 


The coma was found to be highly structured on all scales (jets, 
shells, ion streamers, etc.) and the gaseous component was analyzed 
in situ by mass spectroscopy. Signals at atomic masses of 1 and from 
12 to ~ 55 amu were detected. H 2 O was conhrmed to represent 85% 
by weight of the gas phase, and the likely presence of large organic 
polymeric molecules was indicated. The dust was analyzed by size 
and composition and there was an unexpectedly high fraction of very 
small grains, down to the sensitivity limit of ~ 10“^® kg. Particles 
rich in metals and silicates were found as expected, but particles 
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rich in H, O, C, and N (“CHONs”) were seen for the hrst time and 
were thought to be related to the smallest grains mentioned above 


flKissel et al.L 119861 ). The integrated mass loss experienced by the 
nucleus at this passage, was ~ 0.5% of the total mass of the nucleus. 
The various predicted plasma effects were conhrmed, including the 
existence of a bow shock, and the adjacent interplanetary medium 
was found to be kine matically and magnetically extremely turbulent 


fiFestou et al.l. l2004l) . 


Since the turn of the new millennium, four more comets have 
been visited by various NASA spacecrafts. 

Comet 19P/Bopellv was visited by the Deep Space 1 in 2001 
(jSoderblom et al.l. l2004aj) and the hrst ever dis k-resolved photomet¬ 
ric study was p erformed on a cometary surface (IBuratti et ahl. 12004 : 

20071) . However, since all images were obtained with a 
single hlter, color information for Borrelly was not available. The 
infrared spectrometer onboard Deep Spac e 1 measured the temper¬ 
ature with a single one-dimensional scan f Soderblom et all 2004b ). 

A few yea rs later in 2004, the Stardust mission hew by Comet 
81P/Wild 2 flBrownleel . 120141) . collecting samples from the comet’s 
coma using an aerogel collector and returning them to Earth via 
a sample-return capsule in 2006. The comet samples suggested 
that high-temperature inner Solar System mat erial formed and wa s 


subsequently transferred to the Kuiper Belt f Matzel et ah . 2010[) . 


Moreover, in 2 009, glycine, an am ino acid, was detected in the gath¬ 


ered material f Elsila et ah . 2009[) . 


Comet 9P/ Tempel 1 was v i sited in 2005 and 2011 respectively by 
Deep Impact f[A%earn et ah . 2005) an d Stardust extended mission 
f Veverka et ah 7 120131: Li et al.l . l2013al ). Deep Impact mission was 
al so extended, and in 2010 visited the very active comet 103P/Hartley 


2 flLi et al.l . l2013bt ICronssin et al.l 120131 ) . Main goal of the mission 


was achieved with the detection of water ice ejected from the interior 
of Tempel 1 as a result of the c ollision of the impactor dropped by 
the probe ( Sunshine et ah . 200?!) . Deep Impact, allowed a direct ob¬ 


servation of the inner coma and thu s of the parent spe cies, directly 
linked to the nucleus composition ( Feaga et al.l . 2007 ). Moreover, 
Deep I mpact permitted a hr st map of the surface temperature of a 
comet (jCroussin et al.l . 120071 ). and the h rst spectral detection of ex¬ 
posed water ice on the nucleus surface (jSunshine et al.l . 120061 ) . For 
further details see Chapter [5l 
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1.2 Comets: general view and open questions 


Comets are the most primitive observable objects remaining from 
the era of planetary formation. As such, they provide information 
on the thermophysical conditions of the protoplanetary disk and 
on the formation mechanism for the icy planetesimals from which 
the cores of the outer planets were built. Furthermore, the physi¬ 
cal evolution of cometary nuclei over the past 4.6 Gy must be ex¬ 
plained within the context of any unihed theory of the solar system, 
and comparative studies of cometary nuclei and dynamically related 
bodies should provide insights into the physical and collisional his¬ 
tories of these objects. Through impacts over the age of the solar 
system, cometary nuclei have signihcantly affected the formation 
and evolution of planetary atmospheres and have provided an im¬ 
portant source of volatiles, including water and organic material, to 
the terrestrial planets. Another important motivation for studying 
cometary nuclei is that their bulk properties may dictate what steps 
should b e taken for hazard mitigation in the e vent of a predicted 


collision flLamv et al.l. |2004 iFestou et al.l. 120041) 


Although the long history of cometary science we have still not 
or few information of basic physical parameters like distribution of 
size and mass of the cometary nuclei. Because of they are small 
bodies enveloped in their coma it is possible to have resolved im¬ 
ages only with the most powerful telescopes on Earth, or with is 
situ missions. The principal problem with the ground observations 
is that estimates of the size distribution from optical data require an 
assumed albedo. Many of the data, furthermore, consist of only sin¬ 
gle observations rather than complete rotational lightcurves, which 
adds scatter to the distribution. The cumulative size distribution 
(CSD) of the ECs obeys a single power law with an exponent 1.9±0.3 
down to a radius of 1.6 km. Below this value there is an appar¬ 
ent dehcie ncy of nuclei , poss ibly owing to observational bias and/or 
mass loss f Lamv et ah . 20041) . Regarding the mass we must empha¬ 
size that we still do not have a single measured mass for a cometary 
nucleus, and thus not a single, measured density. A lower limit of 
~ 0.6 g/cm^ for the nucleus bulk density could be found by combin¬ 
ing rotational periods (from 5 to 70 h) and shape dat a when avail¬ 


able to ensure stability against centrifugal disruption fiLamv et ah 
2004 . 
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Determination of the interior structure is limited to what have 
been learned from the Deep Impact mission until the results of 
Rosetta rendezvous and soft lander mission (see section fOl) . Deep 
Impact have been excavate a large crater with a probe that im¬ 
pacted with the nucleus of Comet 9P/Tempel 1 in order to study 
the outermost tens of meters. The expansion rate of the ejecta 
plume imaged during the look-back phase of observations leads to 
an estimate of comet mass (2.3 — 12.0 x 10^^ kg) and a bulk den¬ 
sity of 0.2 — 1.0 g/cm? with large error due to uncertainties in 
the magnitude of coma gas pressure effects on the ejecta parti¬ 
cles in flight. Medium-to-high porosity of the material was also 
establis hed from ejecta plume behavior and from the size of the 
crater ( Richardson et all 2007h . Moreover, implications for inter¬ 
nal stratigraphy of de volatilized material an d water ice of Tempel 1 
were also determined fjSunshine et al.l . 120071) . 


Our knowledge of composition is limited almost entirely to the 
coma, the spectra of nuclei being almost featureless, with the excep¬ 
tion of the water ice features seen on the surface of Tempel 1 and 
Hartley 2. Other than that, the surfaces of comets are known only to 
be very dark, presumably from a combination of particle shadowing 
due to porosity and the inclusion of very dark, carbonaceous mate¬ 
rial as one of the abundant components at the surface. We are thus 
faced with the problem of deciding the extent to which the composi- 
tion in the com a is representative of the composition in the nucleus 
(lA’Hearnl . l2004j j. Interpretation of the chemistry of the coma is fur¬ 
ther limited by the fact that most of the species observed, including 
virtually all the easily observed species at optical and ultraviolet 
wavelengths, are clearly fragments of larger molecules that existed 
in the nucleus. Furthermore, the relative abundances of species vary 
dramatically with heliocentric distance even in a single comet. A 
few cases of variation can be explained in terms of processes in the 
coma, but most are not explained at all and in virtually no case is 
there consensus that we can correct for the variation with heliocen¬ 
tric distance adequately to say something dehnitive about nuclear 
abundances (lA’Hearnl . 120041 ). Our knowledge of the composition of 
the gaseous coma is quite extensive, coming from remote sensing 
at wavelengths from the X-ray to the radio and from in situ mea¬ 
surements. Water is the most abundant constituent of cometary 
ices and its production rate is used for quantifying cometary ac- 
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tivity and for abundance determinations. Other molecules whose 
production rate have been observed to overcome 1% relative to wa¬ 
ter ice are: C 0 ,C 02 ,CH 4 ,CH 30 H,H 2 C 0 ,NH 3 ,H 2 S. Other sul- 
fur, nitrogen, and CHO - bear ing molecules are frequently observed 


f Bockelee-Morvan et al.l . 20041 ) 


The refractory species are much less well known. Remote sensing 
has brought us primarily the silicate feature, including identihcation 
of crystalline olivine, and pyroxene. In situ measurements brought 
us CHON particles, whose presence hlled a major gap in our under¬ 
standing of the overall abundances, since combining these particles 
with the volatiles leads to more or less solar abundances (lA’Hearn . 
2004il . 


Comets contain some of the least-altered material surviving from 
the early solar nebula. Cometary dust may contain both presolar 
particulates and solar nebula condensates. Their structure and min¬ 
eralogy may hold important clues about the chemical and physical 
processes in the early solar system. The key questions in this area 
are (1) whether interstellar ices survived the accretion shock and 
were incorporated directly into comets, (2) whether any chemical 
reactions (either in the gas phase or on grain surfaces) were impor¬ 
tant in that part of the accretion disk in which comets formed, and 
therefore, (3) whether or not the details of the abundances of ices 
in come ts are good con straints on the conditions in the early solar 
system f A’Hearnl . l2004 h 

Finally, among the main issues, there is the role of comets at the 
terrestrial planets. Did comets deliver most of Earth’s water and 
organics? The difference in D/H ratios between comets and terres¬ 
trial ocean water has provided clues in this direction. Measurements 
of this ratio in several Oort cloud comets resulted in a value twice 
as high as that in the Earth oceans, leading to the generally ac¬ 
cepted conclusion that comets are unlikely to be the primary source 
of ocean water. But measurement of the D/H ratio o f 103P/Hartley 
2 resulted to be consistent with terrestrial value (IHartogh et ah 


201111 . 


1.3 Rosetta space mission 

The International Rosetta Mission, approved in November 1993 by 
the Science Programme Committee of the European Space Agency 
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(ESA), is a logical follow-up mission of ESA’s very su ccessful first 
mission, Giotto, to comet IP/Halley (iReinfiardl. 119861) . The prime 
scientific objectives of the Rosetta mission are to investigate the 
origin of our solar system by studying the origin of comets (see 
section It consists of two mission elements, the Rosetta orbiter 
(see Table 11.11 and Figure ll.dh and the Rosetta lander Philae (see 
Table 11.21 and Figure II.Sp . The mission is named after a plate of 
volcanic basalt currently in the British Museum in London, called 
the Rosetta Stone. This plate was the key to unravel the civilization 
of ancient Egypt. Just as the Rosetta Stone provided the key to an 
ancient civilization, the scientihc instruments onboard the Rosetta 
spacecraft are designed to unlock the m ysteries of the olde s t buil ding 
blocks of our solar system, the comets f Glassmeier et al.l . 2007). 

Rosetta will be the hrst spacecraft to orbit a cometary nucleus 
and the hrst spacecraft to hy alongside a comet as it heads towards 
the inner Solar System. This allows for an unprecedented oppor¬ 
tunity to examine from close proximity how a comet, coming from 
deep space, is transformed while heated by the Sun. The lander 
Philae will perform the hrst controlled touchdown of a robotic lan¬ 
der on a comet nucleus, and therefore the mission will provide the 
hrst images taken from the cometary surface and the hrst in situ 
analysis results of cometary material. The Rosetta mission also 
comprises the hrst hybys of two main-belt asteroids. Rosetta is also 
the hrst spacecraft to hy close to the Jovian orbit using solar cells 
as its main power source. This is due to its large solar panels and 
the use of s pecially optimized l ow intensity low temperature (LILT) 


technology (IStrobl et al.l. Il995h . 


The idea that became Rosetta was conceived in the early 1980’s 
even before ESA’s Giotto mission hew by comet IP/Halley, return¬ 
ing the hrst detailed picture of a comet’s nucleus ever seen. The suc¬ 
cess of Giotto meant that plans for a follow-on mission were enthusi¬ 
astically considered. Almost 20 years later, Rosetta was built, tested 
and ready to launch. But disaster struck just one month before the 
planned liftoh. In December 2002, an Ariane 5, similar to the one 
designated to launch Rosetta, failed while lifting a communications 
satellite into orbit. The difficult decision was made to postpone 
the attempt until the launcher failure was understood. This robbed 
the mission of its original target, comet 46P/Wirtanen. While the 
engineers worked to understand and prevent the loss of another Ar- 
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Table 1.1: The Rosetta Orbiter (11 science instrument packages) 


ALICE 

Ultraviolet Imaging Spectrometer 

CONSERT 

Comet Nucleus Sounding experiment 

COSIMA 

Cometary Secondary Ion Mass Analyser 

GIADA 

Grain Impact Analyser and Dust Accumulator 

MIDAS 

Micro-Imaging Analysis System 

MIRO 

Microwave Instrument for the Rosetta Orbiter 

OSIRIS 

Rosetta Orbiter Imaging System 

ROSINA 

Rosetta Orbiter Spectrometer for Ion and Neutral Analysis 

RPC 

Rosetta Plasma Consortium 

RSI 

Radio Science Investigation 

VIRTIS 

Visible and Infrared Thermal Imaging Spectrometer 


floaw* nfM5 



Figure 1.4: The Rosetta orbiter and the scientific instruments onboard. 


iane 5, scientists and engineers searched for a replacement target. 
Eventually they settled on comet 67P/Churyumov-Gerasimenko, a 
somewhat more massive comet than Wirtanen. This led to the 
strengthening of the legs on Rosetta’s lander Philae, to cope with 
the slightly faster landing speed now expected. The launch took 
place on 2 March 2004. 

Rosetta could not head straight for the comet. Instead it began 
a series of looping orbits around the Sun that brought it back for 
three Earth fly-bys and one Mars fly-by. Each time, the spacecraft 
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Table 1.2: The Rosetta Lander (10 science instrument packages) 


APXS 

Alpha Proton X-ray Spectrometer 

CIVA and ROLIS 

Rosetta Lander Imaging Systems 

CONSERT 

Comet Nucleus Sounding experiment 

COSAC 

Cometary Sampling and Composition experiment 

MODULUS 

PTOLEMY Evolved Gas Analyser 

MUPUS 

Multi-Purpose Sensor for Surface 
and Subsurface Science 

ROMAP 

Rosetta lander Magnetometer and Plasma Monitor 

SD2 

Sample and Distribution Device 

SESAME 

Surface Electric Sounding 
and Acoustic Monitoring Experiment 

RSI 

Radio Science Investigation 



Figure 1.5: The Rosetta lander and the scientific instruments onboard. 


changed its velocity and trajectory as it extracted energy from the 
gravitational field of Earth or Mars. During these planetary fly¬ 
bys, the science teams checked out their instruments and, in some 
cases, took the opportunity to carry out science observations coordi¬ 
nated with other ESA spacecraft such as Mars Express, ENVISAT 
and Cluster {http://sci.esa.int/rosetta/). From 8 June 2011 to 20 
January 2014, because too far from the Sun, there would not be 
enough power to keep all the spacecraft’s systems operating. For 
this reason the spacecra ft has been placed in h ibernation, when it 
was in deep space cruise f Glassmeier et ah l . 2007 ). Between wake-up 
and rendezvous with the comet, all twenty-one instruments had to 
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be brought online and checked out. Software was updated and the 
spacecraft had to perform a series of 10 manoeuvres to reduce its 
speed sufficiently to rendezvous with the comet rather than fly by 
it {http://sci.esa.mt/rosetta/). 

The Rosetta mission is comprised of twenty-hve experiments, 
making it unprecedented in scale. The measuremen t goals of the 
Rosetta mission include (ISchwehm and Schnlzl . Il999l ): 


• a global characterization of the nucleus, 

• the determination of its dynamic properties, 

• the surface morphology and composition, 

• the determination of chemical, mineralogical and isotopic com¬ 
positions of volatiles and refractories in the cometary nucleus, 

• the determination of the physical properties and interrelation 
of volatiles and refractories in the cometary nucleus, 

• studies of the development of cometary activity and the pro¬ 
cesses in the surface layer of the nucleus and inner coma, that 
is dust/gas interaction, 

• studies of the evolution of the interaction region of the solar 
wind and the outgassing comet during perihelion approach. 


To achieve these goals Rosetta combines two strategies for char¬ 
acterizing the properties of a cometary nucleus. On one hand, the 
comets evolution along the orbit with decreasing heliocentric dis¬ 
tance will be investigated with the orbiter instruments by monitor¬ 
ing the physical and chemical properties of the nucleus and in situ 
analysis of the near-nucleus environment. On the other hand, the 
Rosetta lander Philae wil l provide ground truth by analyzing the 
nucleus material directly f Glassmeier et al. l . 2007 1. 

The Main ROSETTA Target Comet 67P/Churyumov-Gerasimenko 
is a sho rt-periodic comet f rom the class of the Jupiter Family comets 
(JFC) f Lamv et ah . 2007|l . It was discovered on September 9, 1969 
by Klim Ivanovic Churyumov and Svetlana Ivanovna Gerasimenko 
at the Alma Ata Observatory, Tadchik Republic. JFCs are believed 
to originate from the Kuiper Belt. Since the size distribution in the 
Kuiper Belt is collision-dominated for bodies smaller than about 50 
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Figure 1.6: Steins as imaged by the camera OSIRIS onboard Rosetta during the 
flyby on 5 September 2008, at a distance of 800 km (http://sci.esa.int/rosetta/). 


km, JFCs may actually represent the low-mass fragments from col¬ 
lision events in the belt. They were ejected from the Kuiper Belt 
region by Neptune. Repeated gravitational scattering at the outer 
major planets let them cascade inward. Backward calculations have 
revealed several encounters with Jupiter over the past 200 years with 
the one on February 4, 1 959 being the closest one at just 0.052 AU 
distance from the planet f Belvaev et al.l. 1986 ). This close encounter 
has signihcantly changed the orbit and the perihelion distance from 
2.7 AU before to 1.2 AU after, bringing the comet within reach for 


the Rosetta mission fiGlassmeier et al.l. 120071 ). 


During cruise to its main target the Rosetta spacecraft has also 
visited two asteroids, namely 2867 Steins and 21 Lutetia. Both 
objects are main belt asteroids. 

Steins had been classihed as a rare E-type asteroid based on its 
visual and near-infrared spectrum, and its high albedo. The obser¬ 
vations by OSIRIS fFigureand VIRTIS on Rosetta brought new 
information that could not have been gained from the ground. The 
dimensions of Steins were found to be 6.67km x 5.81km x 4.47km. 
Steins was probably part of a larger differentiated object that had 
broken up. It was later struck by other objects, creating impact 
craters. However, the interior is thought to be a rubble pile and 
the asteroid will eventually break up {http://sci.esa.int/rosetta/). 
Steins probably gained its diamond shape from the YORP effect - 
the Yarkovsky - O’Keefe - Radzievskii - Paddack effect - in which 
photons from the Sun are re-radiated as infrared emission taking 
momentum from the body and altering the rotation rate. In this 
case, the change in rotation rate resulted in some material moving 
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Figure 1.7: Lutetia as imaged by the camera OSIRIS onboard Rosetta on 10 
July 2010. Rosetta flew past at a greater distance, 3162 km, in order to al¬ 
low the full asteroid to appear in the field of view of the scientific cameras 
{http://sci. esa. int/rosetta/). 


towards the equator of asteroid (IKeller et al.l.l2ninh . Surface reshap¬ 
ing could be also due to the YORP effect, accounting for the lack of 
smaller craters on Steins. The effect would have caused landslides to 
hll in the smaller craters. This is the first time that the YORP effect 
has been seen in a main-belt asteroid {http://sci.esa.int/rosetta/). 

On 10 July 2010, ESA’s Rosetta spacecraft flew past asteroid (21) 
Lutetia, one of the largest objects orbiting within the main asteroid 
belt between Mars and Jupiter. 

Discovered in 1852, Lutetia was among the first objects to be 
classified as an M-type (metallic) asteroid, but radar observations 
revealed an unusually low albedo, or reflectivity, that was inconsis¬ 
tent with a metallic surface. Meanwhile, spectra obtained at visible 
and infrared wavelengths found similarities with meteorites known 
as enstatite chondrites and with carbonaceous chondrite meteorites, 
typically associated with C-type asteroids. The flyby of Lutetia 
took place on 10 July 2010, when Rosetta flew past the asteroid at 
a distance of 3170 km. As Rosetta closed on the rotating asteroid, 
the OSIRIS imaging system returned 462 pictures of the illuminated 
northern hemisphere, showing more than 50 per cent of the aster¬ 
oid’s surface. OSIRIS revealed an irregular, cratered object (Figure 
II.7p whose overall dimensions were 121/cm x 101/cm x 75/cm. Its di- 
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verse surface displayed a variety of small, localised features, as well 
as a number of geological features that were on a far larger scale 
than any seen on other asteroids. More than 350 craters were iden- 
tihed with diameters between 600 metres and 55 km and depths 
of up to 10 km. However, the crater density varied considerably, 
the regions with the largest number of craters were clearly older 
than those with fewer craters {http://sci.esa.int/rosetta/). Its an¬ 
cient surface age (determined from crater counting) coupled with its 
complex geology and high d ensity suggest that Lutetia is most likely 
a primordial planetesimal ( Sierks et al.l. 2011 ). This was also con- 
hrmed by t he spectroscopic obser vations performed by the VIRTIS 
instrument f Coradini et ahl . 201 ll ). which have shown no absorption 
features, of either silicates or hydrated minerals. 










Chapter 2 

VIRTIS instrument 
performance analysis 


2.1 VIRTIS 


Spectrophotometry is a very powerful diagnostic tool in remote sens¬ 
ing to study the composition and the physical properties of the sur¬ 
faces of objects under investigation. The amount of solar radiation, 
as a function of the wavelength, scattered from a surface towards 
the observer, is a function of several parameters such as the compo¬ 
sition of the materials making up the surface, their grain size, the 
porosity, the surface roughness and the scattering properties of the 
regolith. On the other hand the thermal emission, is a function of 
temperature and emissivity, which give clues on the composition, 
thermal inertia and active spots on the surface. 

The identihcation of the nature of the main constituents of the 
comets is a primary goal of the R osetta mission, and particularly 
for VIRTIS f Coradini et ah . 200711 (Visual InfraRed and Thermal 
Imaging Spectrometer) that is part of the scientihc payload of the 
Rosetta Orbiter. 

The VIRTIS experiment has been one of the most successful ex¬ 
periments built in Europe for Planetary Exploration. It was devel¬ 
oped under the supervision of Angioletta Coradini. 

The VIRTIS scientihc and technical teams have taken advan¬ 
tage of their previous experience in the design and development of 
spectrometers for space applications. In fact the various groups 
contributing to the VIRTIS experiment from Italy France and Ger- 
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many have been deeply involved in the CASSINI mission with the 
experiments VIMS and CIRS. 

The instrument is composed by a high spectral resolution chan¬ 
nel (-H) and two high spatial resolution channels devoted to spectral 
mapping (-M), covering the spectral range 0.22 - 1.04 fim (Visible) 
and 0.95 - 5.06 pm (Infrared), both with 432 spectral bands. VIR¬ 
TIS will detect and characterize the evolution of specihc spectral 
features, such as the typical spectral bands of minerals and ices, 
arising on the nucleus surface and in the coma. It will also study 
the surface thermal evolution during comet’s approach to the Sun. 
Moreover, it will perform photometric analysis of the surface thanks 
to observations at different phase angles with resolution up to the 
few-meter scale along the mission duration. This work will focus on 
the M channel. Here we intend to show the performance analysis in 
order to characterize the Signal to Noise ratio. 


VIRTIS-M shares the same optical system to analyze two spec¬ 
tral channels (VIS and IR). The angular resolution is 250 prad for 
a single pixel with a FOV of 64 mrad for the whole slit (256 pixels). 
The spectral resolution is 1.8 nm/band for VIS and 9.8 nm/band 
for IR channel. The instrument is designed as a rigid structure in 
order to resist to vibrations caused by the launch, and to maintain 
a high alignment of the two optical heads in flight conditions. Two 
cryocoolers (used for active cooling of the focal planes of the IR 
channels M and H), the boxes of electronics and eight pillars in tita¬ 
nium for the upper part of the instrument are mounted on the pallet 
baseplate (349 x 409 mm). This is located at the base of the struc¬ 
ture and acts as a thermo-mechanical interface with the spacecraft. 
The optical heads are positioned on the upper part of the instru¬ 
ment, in contact with a radiator that keeps the entire instrument in 
radiative equilibrium; the outside of the whole structure is covered 
by a multi-layer insulation. The focus of the telescope, is placed 
at the entrance slit of the instrument, the components of which are 
listed in Figure [2T] ( Filacchionel. 12006 1. VIRTIS-M is equipped with 
a Offner spectrometer, which does not require collimators, lenses or 
beam splitter, and this allows for a substantial reduction of the size 
of the optics (see Figure [2^ . 


The VIS focal plane is constituted by a Thomson CCD, whose 
sensor is cooled to temperatures in the range 150-180 K. The IR 
focal plane is composed by a Raytheon two-dimensional sensor made 
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Figure 2.1: Scheme of VIRTIS: 1. VIRTIS-M; 2. Aperture VIRTIS-M; 3. Cry- 
ocooler VIRTIS-M; 4. VIRTIS-H; 5. Aperture VIRTIS-H; 6.Cryocooler VIRTIS- 
H; 7. Radiator; 8. Isostatic mount in titanium; 9. Pallets baseplate (Selex ES). 



Figure 2.2: Offner spectrometer of VIRTIS-M: I) mirror; 2) diffraction grating; 
3) the entrance slit and electromechanical shutter; 4) infrared focal plane; 5) 
visible focal plane; 6) position of the mirror on visible focal plane (not shown); 
7) mechanical structure (Selex ES). 
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up of photosensitive material (Mercury, Cadmium and Tellurium, 
HgCdTe), with the sensor actively cooled to temperatures 86 - 88 K. 
It is divided into six spectral zones, corresponding to different hlters. 
Since the two focal planes have different thermal requirements, their 
assembly has been realized on orthogonal planes, although very close 
one each other. The CCD receives the diffracted light from the 
grating by means of a mirror that deflects only visible light, while 
the infrared sensor receives the light in a direct way. 

VIRTIS has an internal calibration mode, in order to monitor 
both the spectral and radiometric performance dnring the mission. 
In this mode acquisitions are performed following a sequence with 
different conhgurations of the mechanical components (shntter, cover) 
and two calibration sources, one optimized for the VIS channel and 
one for the IR channel. This provides information on some parame¬ 
ters that can affect the experimental data. In particnlar, the acqui¬ 
sitions with closed shutter permit to obtain the signal of the dark 
current added to the thermal backgronnd, the latter being relevant 
for the IR channel. 


The stored data are saved in a format called “cube”: this name is 
used as VIRTIS-M is in fact capable of recording images at various 
wavelengths. Two dimensions are spatial and the third is spec¬ 
tral. Thus, each cube gives as many mon ochromatic ima.ges as spec¬ 
tral bands. For fu rther details refer to fjCoradini et al.l. 1200711 and 
f Filacchione . 2006). 


2.2 Dark current 

In order to model the dark current in terms of integration time 
{it) and temperature of the focal plane (Tfp) for both VIS and IR 
channel, we assume the following relation: 


Dark{X)DN = a(A) -|- b{\) x it + c(A) x Tpp (2.1) 

The calculation of the coefficients a, b, c, together with the read¬ 
out noise analysis, is performed for the two spectral channels as 
follows. The calculation is performed along the 432 bands (A). The 
dependence on the samples (256 along the slit) is implicit. 
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2.2.1 IR channel 

The acquisitions of dark current investigated in this section have 
all been carried out by VIRTIS during the cruise phase and until 
June 2010 when the S/C entered its hibernation phase. They also 
include dark frames acquired during internal calibration mode. The 
dark frames analyzed include the electronic offset iFigure [273112.4p . 



Bands 


Figure 2.3: Example of acquisition with closed shutter with IR channel (it = 
0.4,Tpp = 87.5K}. 432 bands (horizontal axis) x 256 samples (vertical axis). 
It includes the electronic offset, the dark current and the thermal background. 
The background thermal emission of the spectrometer (T r; 14011) is visible at 
longer wavelengths. The spots are defective pixels. 


During the measurements the IR focal plane is stabilized in tem¬ 
perature by an active cryocooler, while the spectrometer tempera¬ 
ture is regulated by a passive radiator and could vary according to 
the S/C orientation fjCoradini et al.l . 120071 ). 

In order to separate the contribution of the thermal background 
from the dark current, only the hrst 250 bands (1.0 — 3.6 fim) are 
taken into account. 

The dependence of the dark current on integration time and on 
the focal plane temperature is calculated pixel by pixel. The general 
trend of these relations is represented in Fig. 12.51 where the average 
of the pixels of each frame is reported (256 sample x 250 bands). 

While there is a clear correlation between the dark current and 
the integration time, we have observed a negligible dependence of the 
dark current on the focal plane temperature when Tpp = 85 — 88K. 
Many points in the plot are overcoming the 88 K. They correspond 
to acquisitions in internal calibration mode with an integration time 
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Figure 2.4: Equivalent to figure 1^751 for one single sample. The rise above 4/rm is 
due to the thermal background. The different electronic offset along the bands 
caused by the multiplexer readout circuit can also be noted. 
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Figure 2.5: Left panel: Mean Dark current versus integration time. Each point 
is the average of all the pixels in the first 250 bands. The solid line is a linear best 
fit. The Pearson correlation coefficient is used for the analysis. The intercept of 
the line represents the mean electronic offset in DN while the slope corresponds 
to the dark current rate in DN/s. Right panel: same points of left panel as a 
function of the Focal Plane Temperature. 


of 0.5 s. The pattern they show is clearly exponential, althongh 
still negligible if compared to the effect of the integration time. The 
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higher temperature of these points is due to the short time left avail¬ 
able to the cryocooler to cool the focal plane. For larger integration 
times and for temperatures > 88 K this exponential behavior could 
be signihcant. However, the model of dark current does not take 
into account these higher temperatures because the planned acqui¬ 
sitions are supposed to operate within the range where the observed 
correlation is negligible. 

2.2.2 VIS channel 

As in the case of the IR channel also for the VIS channel the acqui¬ 
sitions of dark current taken into account are all those carried out 
by VIRTIS during the overall cruise phase. 



Bands 


Figure 2.6: This dark frame of VIS channel clearly shows an external source of 
light at longer wavelengths. Frames like this are not taken into account in the 
analysis. 

For VIS channel a selection of the dark frames was needed be¬ 
cause in few cases dark current acquisitions are affected by sun light 
contamination passing through a narrow aperture between the ra¬ 
diator and the spectrometer outer chassis. This effect happen only 
for some specihc orientation of the spacecraft. Figure [2T] shows one 
of the dark acquisitions to be discarded for this analysis. 

The amount of dark current is analyzed as a function of the 
integration time and of the temperature of the focal plane. 

Unlike the IR channel the amount of dark current is more affected 
by the temperature of the focal plane than the integration time in 
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the range of possible variation of this quantities (see Figure l2T|) . 
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Figure 2.7: Left panel : Mean Dark current versus integration time. The solid 
line is a linear best fit. The Pearson correlation coefficient is used for the 
analysis. The intercept of the line represents the mean electronic offset in DN 
while the slope corresponds to the dark current rate in DN/s. Right panel: same 
points of left panel as a function of the Focal Plane Temperature. The negative 
correlation observed is due to the low temperature of the sensor in which the 
transfer of the charges is b l ocked by the freezing of the crystal lattice of the 
semiconductor ( Filacchione . 20061) . 


2.3 Readout noise 

2.3.1 IR channel 

The analysis of the readout noise is performed by investigating the 
acquisitions of electronic offset during the internal calibration mode. 
Variations that are not caused by the readout noise have to be 
avoided in the analysis. To this aim the electronic offset is stud¬ 
ied as a function of the temperature of the focal plane, dividing odd 
and even bands. The results are shown in hgure 12.81 and 12.91 

While even bands do not present any signihcant correlation, odd 
bands seem to be correlated with the temperature of the focal plane. 
In order to evaluate the readout noise, the contribution of the tem¬ 
perature has to be eliminated choosing frames acquired at the same 
temperature. Many (50) acquisitions of the offset with a focal plane 
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Figure 2.8: Electronic offset for even bands (IR channel) in function of the 
temperature of the focal plane. There is not significant correlation. 



Figure 2.9: Electronic offset for odd bands (IR channel) in function of the 
temperature of the focal plane. There is a negative correlation. 

temperature of ~ 87 K (the most frequent) are selected for the anal¬ 
ysis. Then the standard deviation is performed pixel by pixel. The 
average of the standard deviations over all pixels is a measurement 
of the readout noise in digital number. The result for IR channel 
is 2.6 DN, equivalent to 175 photoelectrons. It is assumed to be 
independent from instrumental parameters. 





46 CHAPTER 2. VIRTIS INSTRUMENT PERFORMANCE ANALYSIS 

2.3.2 VIS channel 

Here the same procedure of section 12.3.11 is followed. For the visible 
channel there is no difference for the electronic offset between odd 
and even bands. The correlation of the electronic offset with the 
temperature of the focal plane is negative as for the dark current (see 
section [2. 2. 2p . In order to analyze the readout noise 50 acquisitions 
with the same focal plane temperature (~ 155 K) are analyzed. The 
resulting readout noise is 1 DN, equivalent to 65 photoelectrons. 



Figure 2.10: Mean offset of VIS channel vs Temperature of the focal plane. The 
correlation is negative. 


2.4 Thermal background 

The IR channel is further affected by the thermal background which 
has to be taken into account for the calculation of the noise and the 
saturation limit. 

We model the thermal background as a function of the integration 
time {it) and the temperature of the spectrometer (Tg): 


Background{X)DN = K{X) x B{\,Ts) x it x ITF{X) (2.2) 


where: 
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B{X,Ts) is the internal thermal radiance depending on wave¬ 
length and spectrometer temperature; 

K{X) accounts for the hlter transmittance, the effective area 
of the spectrometer seen by the detector and other possible 
contributions. 


ITF(X) is the instrument transfer function in {DNm?umsterad/W /s 


flFilacchionel.l2006l:lFilacchione et al.l. 120061: iMigliorini et al.l.l2013l ). 


The characterization of the K{X) parameter is crucial for mod¬ 
eling the thermal background. K{X) can be derived from measured 
acquisitions of thermal background with known Ts and it by invert¬ 
ing the Eq. 12.21 Acquisitions with different Ts and it are taken into 
account to perform this analysis. Figure [2.111 shows a simulation of 
the Thermal Background with different spectrometer temperatures. 



Wavelength (|im) 


Figure 2.11: The simulated DN as a result of four different temperatures of the 
spectrometer, with an integration time of 1 s. Thermal Background affects the 
spectrum at decreasing wavelength for increasing temperature. The pixels are 
saturated where the line is missing. This simulation is performed for the central 
sample of the slit. The result is similar for the other samples. 


2.5 Signal to noise ratio modeling 

The S/N Simulator is a software tool able to calculate the expected 
instrument signal to noise ratio (S/N) for different input signals and 
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observing conditions. Only Poissonian noise is taken into account 
during the calculation, as other kind of noises are assumed negligible 
for our purpose. 

The calculation performed to obtain the S/N depends on the kind 
of signal to treat: measured or simulated. 

2.5.1 Measured signal 

A correct characterization of the noise for measured signals is im¬ 
portant in the retrieval phase of the surface properties as discussed 
in chapter m 

Prior to the calibration, a raw signal is recorded in Digital Num¬ 
ber (DN), directly convertible to PhotoElectrons (PE). From the 
latter the Poissonian noise and the resulting S/N are straightfor¬ 
ward: 


Noise{\)pE = -s/Total signal{X)pE + readout noise‘s (2.3) 

S/N{X) = R{X)/Noise{X) (2.4) 

TotalSignal{X) includes the signal of the target and all other 
sources of signal such as the dark current and thermal back¬ 
ground. It is obtained by subtracting the electronic offset to 
the measured acquisition. 

R{X) is the signal of the target. It is isolated from other sources 
of signal by subtracting an acquisition with closed shutter to 
the acquisition with open shutter, being the former sufficiently 
close in time to the latter. 

Moreover, the tool checks for saturated pixels, if Total Signal + 
electronic offset overcomes the saturation limit of the instrument. 

2.5.2 Simulated signal 

The analysis of the expected noise for a given simulated signal is 
mandatory to plan the observations. An example is presented in 
section 13.41 where the study of the detectability of minor compo¬ 
nents is analyzed as a function of the integration time and other 
parameters which affect the instrumental noise. 
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In order to calculate the S/N, for simulated signals we have to 
simulate the Total Signal as well, which include dark current and 
thermal background. 


Total signal{\) = R{\) + Dark{X) + Background{X) (2.5) 

From Eq. 12.51 we can then compute Eq. 12.31 and 12.41 to derive the 
S/N. 

The terms of Eq. 12.51 depend on the integration time and the 
temperature of different parts of the instrument. In order to perform 
simulation of the dark current and thermal background we have 
analyzed the DN acquired by the instrument in different conditions 
for both Visible and Infrared channels as explained in sections 12.2[ 
12.31 and 12.4[ The simulation of the target signals R{X) is addressed 
in chapter [3l 

2.5.3 Testing the S/N simulator 



Wavelenglh 


Figure 2.12: The simulated S/N (blue line) is over plotted to the S/N extracted 
from the visible lamp signal (black line) acquired during internal calibration 
mode with IR channel. The integration time is 0.1 s. The average overestimation 
of the simulator is 30%. This comparison also represents a satisfactory test for 
the saturation of IR channel (where the S/N is set to 0). The band at ^ 2800 
nm is the limit before which the signal is saturated. This is what emerges from 
the real data and the simulation. 
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Figure 2.13: The simulated S/N (blue line) is over plotted to the S/N extracted 
from the infrared lamp signal (black line) acquired during internal calibration 
mode with IR channel. The average overestimation of the simulator is 40%. 
Integration time = 0.5 s. 



Waveleiiglh 


Figure 2.14: The simulated S/N (blue line) is overplotted to the S/N extracted 
from the Lutetia spectral cube (black line) in IR channel. The integration time 
is 0.7 s. The average overestimation of the simulator is 40%. 


To test the accuracy of the tool in reproducing the S/N, a method 
has been developed to extrapolate the S/N from real acquisitions of 
the VIRTIS instrument, which are to be compared with the results 
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Figure 2.15: The simulated S/N (blue line) is over plotted to the S/N extracted 
from the visible lamp signal (black line) acquired during internal calibration 
mode with VIS channel. The integration time is 1 s. The average overestimation 
of the simulator is 15%. 



Figure 2.16: The simulated S/N (blue line) is over plotted to the S/N extracted 
from the infrared lamp signal (black line) acquired during internal calibration 
mode with VIS channel. The integration time is 20 s. The average overestima¬ 
tion of the simulator is 10%. 

of the simulator. 

The method is aimed to study the signal poissonian fluctuations 
compared to the average, excluding fluctuations due to changes in 
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Figure 2.17: The simulated S/N (blue line) is over plotted to the S/N extracted 
from the Lutetia spectral cube (black line) in VIS channel. The integration 
time is 1 s. The simulator underestimates the S/N in average by a factor 15%. 
This underestimation of S/N is due to the underestimation of the signal. In 
fact the VI S channel is affect ed by contamination from the high orders of the 
IR channel ( Filacchionel l2006ll . and this contamination is not well characterized 
by the ITF. This is a further source of error that is not taken into account in 
this analysis. 


radiance. 

A number N of radiances spectra, from the same source, that 
must remain as stable as possible, have been investigated. The 
mean radiance is the reference against which the fluctuation over 
N radiances is calculated (for each band). The fluctuation due to 
changes in radiance are overcome by adapting all the spectra to the 
mean spectrum by a running box of the size of 20 bands: the signal 
of each band (central band of the box) is divided by the average 
value relative to the box bands and multiplied by the average of the 
box bands of the mean radiance. This procedure is repeated for each 
one of the N radiance spectra. Residual fluctuations of the resulting 
spectra for each band are due to poissonian noise, and from that the 
S/N is calculated. 

As stated above, the source has to be as stable as possible. Three 
sequences of acquisitions have been identifie d that comply wi t h this 


requirement: two from calibration lamps flMelchiorri et al.l. 120031) 


whose signal is acquired during the internal calibration mode, and 
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a sequence of acquisition of Lutetia asteroid, during its flyby. 

The available acquisitions are: 80 for Lutetia’s frames, that is a 
full sequence; 50 for internal lamps, that is a sequence maintaining 
(enough for our purpose) the same focal plane and spectrometer 
temperatures. 

The quality of the result depends on the available statistics. The 
larger is the number of spectra the more reliable will be the result. 
An average of the resulting S/N along the samples could compensate 
the poor statistics. 

The available samples are: 12 for Lutetias frames, that are the 
samples always hlled with asteroid surface during the sequence; 256 
for internal lamps acquired in visible channel, that is the whole slit; 
128 for internal lamps acquired in infrared channel, that is the half 
slit, avoiding the p art affected by a defect of the order sorting hlter 
(FUacchione, 2006 1. 

Figures 12.121 12.131 and 12.141 show the result of the comparison 
for IR channel, while VIS channel is described by hgures 12.15112.161 
and 12.171 As expected, in general the simulator overestimates the 
S/N because it does not take into account other sources of noise 
than the Poissonian noise. To compensate this overestimation, in 
the analysis preformed in section [21 a conservative approach is used, 
lowering the resulting S/N by 40 % for IR channel. 

Looking at the results, for VIS channel a correction is not needed. 
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Chapter 3 

Spectral simulation 


In order to simulate how VIRTIS-M handles the signal coming from 
the nucleus’ surface of the comet, we have to dehne the physical 
properties of the surface and apply a radiative transfer model to 
calculate the spectral radiance entering the instrument. The simu¬ 
lation of the spectra is also dualized to test the developed tools for 
the spectrophotometric analysis discussed in next sections. 


3.1 Hapke model 


The Hapke model flHapkel . Il993l. l2012bl) is used both to simulate the 
spectra and to retrieve quantitative information on the parameters 
describing the surface physical properties. The most recent version 
accounts for the effects of regolith porosity on the overall redectance. 
Recent work s have established its greater reliability than the previ - 
ous version (IHelfenstein and Shepardl . 120111: ICiarniello et al.l . 120141) . 

In Hapke’s theory the surface bidirectional redectance can be 
expressed as: 


r {i, e,g)=K{oj/ dvr) [/xoe/(/ie + l^oe)] 

{p{g)[l + BsoBsig)] + H{fioe/K)H{fijK) - 1} (3.1) 
[1 Bco Bcig)] S{i,e,g,e) 

where: 


Poe, Pe- edective incidence and emission angle cosine. They are 
the cosines of the angles corr ected by the ef fect of the surface 
roughness (see section 12.4 in Hapke f 2012bl) . 
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oj\ single scattering albedo (SSA). It contains the spectral in¬ 
formation. Once the optical constants and particle diameter 
are hxed it is possible, following Hapke’s model, to compute 
single scattering albedo for a given type of particles. It can 
assume values in the range 0-1. 


p{g)'. single scattering phase function {g is the phase angle) 
modeled with Heyney-Greenstein function with 2 parameters 
b,v: the hrst one describes the angular width of both forward 
and back scattering lobes (0-1 r ange), while the second one de¬ 


scribes their relative amplitude flHenvev and Greensteinl . 11941 
Domingue and Verbiscer . lOOTt) : 


p(^) = - (1 —-[. 

2 (1 - 2bcos{g) + 62)3/2 ^ 

2 (1-1- 2bcos{g) + 62)3/2 


(3.2) 


it must be noted that p{g) cannot be negative. Thus in depen¬ 
dence of the 6 values, y has p hysical meaning within the range 
permitted by ( Hapkel . l2012af) : 


n < 


(1 + 362 ) 
6(3 + 62) 


(3,3) 


Bso, Bs{g), Bqo and Bc{g) describe the observed non-linear in¬ 
crease in reflectance towards small phase angles due to two 
mechanisms: the Shadow Hiding Opposition Effect (SHOE) 
and the Goherent Backscattering Opposition Effect (GBOE). 
In principle the amplitude of both {Bso, Bco) should not be 
greater than 1. Bs{g), and Bc{g) are linked to the angular 
amplitude hs and he- To model them we choose the approxi- 
m ated eq u ations (respectively 9.22 pag. 232, and 9.43 pag. 244 


m 


Hapke 


gn at ions (^re 

3 fl2012b[) ). 


K: porosity parameter. It is linked to the hlling factor 0, that 
is the total fractio n of the volume occupied by the particles (Eq 
7.45b, pag. 167 in Hapke f 2012bl )). and it is meaningful up to 
the critical point K ~ 1.9, corresp onding to (t> < 52% where 


coherent effects become important f Hapkel . 2012bl) 
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H: Ambartsum ian-Chandrasekhar function (Eq. 8.56, p.204 in 
Hapke f 2ni2bl )). it is related to the multiple scattering, which 


involves a; in a non-linear way. 

S{i,e, g,9): shadowing function, is the correction due to surface 
roughness. Its value is always less than 1 and decreases with 
increasing roughness parameter 9, which is an averag e slope 
of the fa cets composing the surface (see section 12.4 in iHapke 

fl2012bh h 

We remand to Hapkel (2012b) for further details. 

To obtain the equivalent radiance spectra we need to model the solar 
irradiance and the thermal emission: 


R = J/D^r + B{T,eeff) 

where: 

r : bidirectional reflectance (Eq. 13.ip 

J : solar irradiance measured at 1 AU 
{http://rredc.nrel.gov/solar/spectra/amO/modtran.html) 

D : heliocentric distance in AU 
B : thermal emission 
T : temperature 


(3.4) 


eg// : effective integrated emissivity fiDavidsson et al.l. l2009h 


3.2 Dark terrain definition 

To obtain realistic spectra of the surface we need to define a main 
component which for comets has typically a very low albedo (re¬ 
ferred to hereafter as Dark Terrain). The definition of the SSA for 
the Dark Terrain is crucial for the hnal model of reflectance com¬ 
ing from the mixing with minor components. This in turn dehnes 
the simulated detectability for the VIRTIS instrument, and so the 
conditions (environmental and instrumental) required to perform 
abundance measurement. To dehne the SSA of the Dark Terrain 
we take advantage of the data obtained by the HRI-Deep Impact 
spacecraft, that encountered Comet 9P/Tempel. 
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Ground observations suggest that comet Tempel 1 has very simi¬ 
lar phot ometric prop e rties to that of CG 67/P, target of the Rosetta 


mission f Lamv et ah . 20061 ) 


We analyzed publicly available data (http://pds.jpl.nasa.gov/) 
obtained by the HRI (High Resolution Imager) onboard Deep Im¬ 
pact consisting of a telescope servi ng both a camera and a n imaging 


spectrometer in the infrared range ( Hampton et ah . 20051) . We have 


selected data taken before the impact, where the nucleus is resolved 
and has a maximum spatial resolution of 120 m per pixel. We use the 
data of the spectrometer in t he range 1.2 - 3.5 /i??r which present an 
higher radiometric accnracy (IKlaasen et al.l . l2008l. 120131) . To obtain 
the reflectance from the radiance, we modeled it as the snm of two 
contribntes: the reflected snnlight (1.2 — 2.5 pm) and the thermal 
emission (2.5 — 3.5 pm), setting as free parameters the tempera- 
tnre, the eff ective integrated e missiy ity and the spectral reddening, 
similarly to iDavidsson et al.l (120091) . The data of the camera are 
available for seven different bands from 0.35 to 0.95 pm. These im¬ 
ages are rescaled to the same spatial resolution of the hyper-spectral 
cubes obtained from the spectrometer. In this way we are able to 
extend the spectra of each pixel to the visible wavelengths. This task 
is straightforward because all the data we have used were acquired 
in a short time interval compared with the comet’s rotation period. 
Therefore, the dataset we are analyzing contains only observations 
taken on the same nncleus hemisphere. Since the radiometric accu¬ 
racy of the VIS channel is higher than the IR, we have scaled the 
latter to the former before bridging the spectra of each pixel. Ont of 
statistics pixels are then removed by applying a the despiking hlter 
discussed in section [6.2.11 on the whole dataset. 

The hyperspectral cubes containing both visible and infrared 
data are powerfnl tools to recognize icy regions. After having as¬ 
signed thresholds for the reflectance at visible wavelengths and for 
the 2.0 pm band dep th is possible to ide ntify the ice-rich nnits. 


The resnlts ob t ained ( Raponi et ah . 20131) are similar to that of 


Sunshine et al.l (120061) . To dehne a SSA of the Dark Terrain we 


need to exclude the icy units from the analysis. Taking into acconnt 
the Hapke model (Eq. 13.11) . we can derive a SSA from the average 
spectrnm of reflectance of the non-icy nnits. However, because the 
available observations are all taken at a single phase angle (63°) we 
have to assume the value of the parameters in the model. At this 
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specific phase angle it is assumed that the opposition effect (sig- 
nihcant at small phase angles) and the shadowing caused by the 
roughness (signihcant at high phase angles) are negligible. More¬ 
over, the porosity parameter and the single particle phase function 
are hxed to 1. In this way the SSA has been derived by inverting 
the Eq. 13.11 and a linear interpolation of it is performed (see Fig. 
OandEUD. 



Figure 3.1: The defined SSA of the Dark Terrain is plotted in black. Two 
different mixing are over plotted. The Dark Terrain is mixed with water ice, 
carbon dioxide and methanol, with different abundances. 


3.3 Synthetic spectra 

The various compositional end-members present on the surface can 
be mixed in different ways. Here we selected two mixing modalities: 
areal and intimate mixing. 

In the case of areal mixing, namely the surface covered by patches 
of different composition, the SSA of each depends on the spectral 
properties of a single component. The resulting reflectance is then 
a linear combination of the reflectances of the various components. 
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weighted for the relative end-members abundances. 

In the case of intimate mixing particles of different composition 
are mixed together (“salt and pepper” mixture) and the hnal SSA 
is a linear combination of the SSA of the different components. 

In modeling the comet surface we take into account different mi¬ 
nor components (and the respective optical constants) which could 
be found in icy form on the surface. They are mixed with the Dark 
Terrain as dehned in previous section. 

The SSA of the minor components simulated in this work comes 
out by the following optical constants: 


water ice (H 9 O) flWarrenl . Il984l : iMastrapa et al.l . 120081 12009 


Clark et ah . 2012[) 


carbon dioxide {CO 2 ) I Qnirico and Schmittl . 1997h (GhoSST 
service http://ghosst-prod.obs.ujf-grenoble.fr/) 

methanol {CH 3 OH) ( Trotta and Schniiti] . 199fil ) (GhoSST 
vice http://ghosst-prod.obs.ujf-grenoble.fr/) 


ser- 


To simulate a real spectrum we add noise to the theoretical spec¬ 
trum. 

The base for its calculation is the S/N(A) coming from the S/N 
simulator (section l2.5p . The S/N simulator converts the simulated 
radiance R into digital number (and photoelectrons) as a function 
of the integration time: 


R{X)dn — X ITF(^X) X it (3-5) 

This quantity is used in Eq. 12.51 for the calculation of the S/N. 

The error bars for radiance and reflectance spectra are: 


CTRadiX) = Radiance/{S/N) (3.6) 

(jRefi{X) = Radiance/{S/N)/Solar Irradiance (3.7) 

Therefore, for a given radiance spectrum in input and given in¬ 
strumental and observational conditions the S/N simulator gives as 
output the error bars of the spectrum. They could be considered 
as the standard deviations of the random fluctuations of the signal 
around the value of each band. 
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3.4 Minor component detectability 


The detection of the spectral bands is the hrst step in order to detect 
the composition of the surface, from which follows the retrieval of 
the surface model parameters as the abundance of the end-members 
and the grain size. 

Since the detection of the spectral bands depends on environ¬ 
mental and instrument conditions (in addition to the composition 
of the nucleus), we need to summarize the possible observation sce¬ 
narios, in order to obtain a general view, and to identify the optimal 
integration time as a function of the case considered. 

In modeling comet’s surface we have assumed three different mi¬ 
nor components (and their respective optical constants) which could 
be found in icy form on the surface: water, carbon dioxide, and 
methanol. The presence of the ice is expected in regions cold enough 
to prevent sublimation. In particular in the early morning we could 
observe frost ice formed in night time as a result of the deposition of 
material ejected by active spots. The ice could also be found in re¬ 
gions with higher temperature if it is physically sep arated from non 


icy co mponents, as assumed in the case of Tempel 1 fISunshine et ah 

2nnfih . 


Four different compositions of the nucleus are simulated: areal 
and intimate mixing, both with high and low abundance of the com¬ 
ponents selected (see previous section): 

1. Intimate mixing with 1% of water ice, 1% carbon dioxide, 1% 
methanol, 97% dark terrain 

2. Areal mixing with 1% of water ice, 1% carbon dioxide, 1% 
methanol, 97% dark terrain 

3. Intimate mixing with 5% of water ice, 5% carbon dioxide, 5% 
methanol, 85% dark terrain 


4. Areal mixing with 5% of water ice, 5% carbon dioxide, 
methanol, 85% dark terrain 

In order to simulate typical VIRTIS-M observations of Churyumov- 
Gerasimenko comet from Rosetta’s orbit we have considered the 
following observation scenario: 

• e = 0", because the geometry of view is supposed to be in nadir 
direction. 
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• i = 80°, the output strongly depends on this parameter and 
this choice represents an observation in the early morning (local 
time). 


Grain size = 50 pm. It is an intermediate value with re- 
spect to the results obtained in our Tempel 1 data analysis 
( Raponi et al.l. 2013h . 


Single particle phase function backscattering: Heyney-Greenstein 
function parameters: b = 0.2, v = 0.73, being th eir value linked 
with the “hockey stick relation” { Hapkel . 2012a[l 

Roughness parameter 9 = 16°, a. s resul ting from Tempel 1 pho¬ 
tometric data analysis I Li et ahl . 2006h . 

Porosity and opposition effects assumed negligible. 


• Heliocentric distance in AU = 3.2, 2.5, 1.5. These values are 
related to different phases of Rosetta’s mission. 


The temperature of the nucleus is linked to the heliocentric distance 
with the simple equation representing the energetic equilibrium: 


5(1 - A) p/D^ = (1 - ~^h<ysB (3.8) 

Where: 

S is the solar constant; 

A is the bond albedo, calculated from spectral information of 
the components; 

D is the heliocentric distance; 

eh is the integrated emissivity, calculated from spectral infor¬ 
mation of the components (see section [5.1.2p : 

^ is the self-heating parameter. It is assumed equal to 0.5 

value with respect to those obtained by 
for Tempel 1 surface; 

a SB is the Stefan-Boltzmann constant; 

T is the temperature of the nucleus surface. 


that is an intermediate 


Davidsson et al.l 020091 1 
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This is an approximated expression because it does not consider 
the contribution of the thermal inertia of the comet’s nucleus, but 
for our purpose it is sufficient in order to provide different plausible 
temperatures of the nucleus’ surface. 

For a given heliocentric distance, and a radiance spectrum in in¬ 
put we have simulated the reflectance spectra and the error bars 
associated to each band, being the latter resulting from the calcu¬ 
lated S/N (Fig. ra . 

In order to determine the achievable accuracy in the characteriza¬ 
tion of spectral features it has been developed a tool that calculates 
the error associated to each band area thanks to the information 
coming from the S/N simulator. The band area is calculated with 
the formula: / 1 — ref lectance/continuum dX. 

The used algorithm performs the calculation taking into account 
possible variations in the shape and wings position of the absorption 
bands due to different mixing methodologies. 

The following diagnostic features are considered : 

• H 2 O ice (1500 nm) 

• CO 2 ice (4200 nm) 

• CH^OH ice (3600 nm) 

The plots in Fig. 13.3113.4113.51 show the error associated to the 
band area as a function of the integration time for low (left) and 
high (right) abundances, and for intimate (up) and areal (down) 
mixing. Each hgure shows a different heliocentric distances. The 
end-members are mixed with the Dark Terrain. 

As expected, factors that increase the detectability of a spectral 
feature for a given end-member are: 1) longer integration time; 2) 
higher end-member abundance; 3) shorter heliocentric distance. On 
the other hand, if the signal is too high there is the possibility to 
reach saturation. Identifying an optimal integration time is therefore 
necessary for balancing the various factors. An integration time is 
considered optimal if it is far from the time at which saturation 
begins and if it permits to correctly detect the end-members. 

In the hrst phase of the mission the comet is far from the Sun. 
The signal of both reflected sunlight and thermal emission is low. 
The error associated to each band area is high for short integration 
times. In general an integration time > Is is needed for the detection 
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of water ice, but a time > 5 s is required to properly detect carbon 
dioxide and methanol. However, close to the sun (1.5 AU) 5 s is 
also a limit beyond which saturation is possible, due to the high 
signal of thermal emission and reflected sunlight, which prevents 
the detection of the bands respectively of carbon dioxide and water 
ice. On the other hand the error on the band area is decreasing 
because of the higher signal from the nucleus surface. 

Although the S/N Simulator is capable to handle variable spec¬ 
trometer’s temperatures, in the study shown in this section we have 
assumed a nominal temperature of 140 K which is in fact the mean 
temperature recorded during the mission. However, during some 
limited phases of the mission, the spectrometer’s temperature has 
risen up to ~ 150 A'. For this reason we have verihed how the results 
are affected by an increment of the spectrometer’s temperature. We 
repeated the same calculation with Ts = 145 K and Ts = 150 K. In 
general the errors on the band area slightly worsen. However, the 
main issue is represented by the saturation of pixels that begins at 
wavelengths gradually smaller with increasing spectrometer’s tem¬ 
perature. In some cases, the CO 2 band is no longer detectable for 
Ts > 145 K and integration time > 5 s, even at 3.2 AU where the 
thermal emission from the nucleus surface is lower. 

As stated in previous section the variation of the focal plane 
temperature in the IR channel produces negligible variations in the 
total signal, at least for the range of temperatures recorded during 
the mission. 
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Figure 3.2: Retrieved radiance (up panel) and reflectance (down panel) with er¬ 
ror bars. In this example the simulated composition is: intimate mixing with 5% 
water ice, 5% carbon dioxide, 5% methanol, 85% dark terrain. The heliocentric 
distance is fixed to 2.5 AU, from which the resulting temperature of the surface 
is 189 K. The integration time is fixed to 2 s. The error bars are calculated as 
explained in section 4.3. The error bars for radiance spectrum are enhanced by 
a factor 10 to make them clearly visible. Other parameters are fixed as stated 
in the text. 
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3.2 AU 



02468 10 02468 10 

Integration Time (s) 


Figure 3.3: Relative error on the retrieved band area versus different integration 
times, for heliocentric distance of 3.2 AU. The four panels represent different 
mixing modalities: intimate (up), areal (down), 1% abundances (left), 5% abun¬ 
dances (right). 







Band Area Error (%) 
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2.5 AU 



02468 10 02468 10 

Integration Time (s) 


Figure 3.4: Same of Figure [3^ for heliocentric distance of 2.5 AU. 
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1.5 AU 



Figure 3.5: Same of Figure l3.31 for heliocentric distance of 1.5 AU. The missing 
points represents saturated pixels in correspondence of the absorption bands of 
the end-members. 






Chapter 4 


Model-Parameters 

extractor 

4.1 Algorithms for spectrophotometric analysis 

In general the observed data need to be compared with a model to 
retrieve the surface’s properties. Referring to the Hapke equations 
fl3.1l and 13.4p we can distinguish three sets of parameters, which 
can be retrieved separately because they concern different physical 
processes: 

• parameters describing the thermal emission (T and e ); 

• parameters describing the position and the shape of the absorp¬ 
tion bands (end-members abundances and grain size) to define 
the SSA; 

• parameters describing the surface physical properties and struc¬ 
ture: 

{p{g),9,K,Bso,Bs{g),Bco,Bc{g)). To retrieve them we need 
several spectra taken at various phase angles. 

The parameters in each set should be retrieved simultaneously to 
avoid hctitious constrains. The retrieval procedure has to search for 
the minimum of the Xr variable, namely the best £t between the 
model (r™) and the observed (r") radiance or reflectance. 

XR = {SAl(r”-rr)/ff.t}/r>OF (4.1) 
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Where a is the calculated noise relative to the observed data. 
DOF are the degrees of freedom, equal to the bands of the spectrum 
involved in the retrieval process, reduced by the numbers of the 
parameters to retrieve. 


The Xr variable is defined in a N-dimensional space, where N 
is the number of the free parameters to retrieve. So basically the 
three sets of parameters define three different spaces. We could need 
many parameters (many dimensions of the space) to fit the data. For 
example many end-members could be required to £t the absorption 
bands of a spectrum, and up to 9 parameters are required to £t 
the phase function with the Hapke model. To perform this task we 
can take advantage of some curve-fitting algorithm generally used to 
solve non-linear lea st square s prob lems, in particular the gradient- 
descent al gorithm (Snvman . 20051) and the Levenberg-Marquardt 
algorithm (iMarkwardt . 20091) . being respectively a first-order and a 
second-order optimization algorithm. 

Like other numeric minimization algorithms, these algorithms 
are iterative procedures. To start a minimization, the user has to 
provide an initial guess for the parameter vector. In cases with 
only one minimum, an uninformed standard guess like (0,0,...,0) 
works hne; in cases with multiple minima, the algorithm converges 
to the local minimum. For this reason more than one first guess are 
adopted, trying to explore all the range of possible values for each 
parameter. The absolute minimum is considered as the best £t if 
the Xr value has a significant lower value than other local minima. 


The best fit gives back a more reliable result if the information 
of the noise coming from the S/N simulator is used in the definition 
of the (Tj in Eq. 14.11 (see Fig. 14.3p . This can be done for every real 
spectrum acquired by the instrument. This happens because each 
band is better weighted during the best fit process. Moreover, the 
resulting Xr function coming from the best fit could be compared 
with the theoretical Xr distribution, whose average is 1. This can 
give a better evaluation on the goodness of the model retrieved. 

The spectral analysis procedure takes as input the observed data, 
the geometry of the shape model, the calculated S/N of the signal, 
and an optical constant database. Information on geometry and 
S/N ratio are not mandatory: the absence of geometry information 
means we can model only the shape and the position of the absorp¬ 
tion bands. In this case the tool can calculate a spectral fit, but a 
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photometric analysis is not possible. In absence of information on 
the S/N the tool can use another reasonable weight (a) in Eq. 14.11 
like the measured signal. Without these information the retrieval, 
even still possible, loses reliability. 



Wavelength (|im) 

Figure 4.1: As a first guess the reflectance spectrum is modeled taking into 
account the dark terrain and the optical constants of the water ice and the 
carbon dioxide. Up: the synthetic spectrum and the resulting best fit (clear 
line). The reduced chi-square results to be 3.6. Down: the residuals spectrum 
smoothed with a 10 bands running box. Only the residuals having the signal 
overestimated by the model are shown to enlighten absorption bands we are 
missing. We consider significant a residual if it overcome 2 a. We can notice a 
significant residual at ^ 3.5 /rm. We can argue we are missing the absorption 
band of methanol. 

A simple example of the spectral analysis procedure, with simu¬ 
lated spectra, is shown in Figure 1441 - 14.21 The specihc composition 
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Wavelength ()J,m) 

Figure 4.2: We repeat the procedure of Figure 14.11 including methanol. The 
resulting reduced chi-square is 1.1. The bottom panel show the residuals spec¬ 
trum keeping the same range of Figure 03] for the sake of clarity. There are no 
significant residuals. We can argue we have achieved a good model. 


and viewing geometry is the same shown in hgnre 13.21 being the 
relative instrnmental, geometric and compositional parameters as 
stated in section 13.41 The respective noise is also simnlated. 

The resnlting reflectance spectrnm is analyzed as discussed. Thanks 
to the analysis of the residuals ([r° — 'r™]/cTi) a final model is ob¬ 
tained. The fitting procedure have given as output the same abun¬ 
dance and the grain size used to simulate the spectra, showing the 
functionality of the method. 

With the availability of an adequate optical constant database, 
and the capability of the fast algorithms to handle many end-member. 
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we can perform an automatic procedure to select the end-members 
as a function of the position of the signihcant residuals. 

Results from the spectral £t allow us to determine the mixture 
and the grain size. Their knowledge enables us to compute and £x 
the single scattering albedo u. Once we have a collection of spec¬ 
tra of the same surface region observed at different phase angles we 
can model the phase function with a best fit with the Hapke model. 
The possibility to £x u reduces the complexity of phase function £t 
procedure and allows the decoupling of spectral effect from photo¬ 
metric ones. A £t of the phase function could be performed at each 
wavelength. However, with the knowledge of u we can perform a 
phase function £t taking into account the whole spectrum, or part of 
it, allowing a better retrieval of the parameters thanks to the larger 
statistics. In this way the £t of the phase function is performed 
after the spectral £t. However, the procedure could be inverted. A 
phase function £t could be performed at each wavelength before the 
spectral fit. In this case the single scattering albedo a; is a free pa¬ 
rameter to retrieve. This approach is followed for Lutetia analysis 
showed in Chapter 16.11 

Once we have a spectrum of the retrieved single scattering albedo 
it can be modeled by a spectral £t, or it can be assumed as end- 
member for spectral unmixing of nearby regions which present sim¬ 
ilar spectral features. 

4.2 Parameter errors evaluation 

4.2.1 Errors from signal uncertainties 

We can be interested to link the error bars of the spectra to the 
error of the retrieved parameters. An intuitive method to estimate 
the total error on the retrieved parameters due to the uncertainties 
on the input data is given here. After the retrieval procedure we 
obtain values of spectrophotometric parameters by which we have a 
model of the spectrum. The model could be used to generate many 
synthetic spectra adding random poissonian noise, thanks to the 
information of the a as estimated by the S/N simulator. For each 
synthetic spectrum generated in this way, the retrieval procedure 
is repeated, giving as a result different values than those obtained 
originally. At the end we have a gaussian-like distribution of values 
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for each parameter. Taking the 95% of the distribution would be 
a conservative estimation for the error bars of the parameters. An 
example is shown in Figure 031 



Retrieved Water Ice Abundance (%) 

Figure 4.3: Four thousand noisy spectra are simulated, and the plots show 
the probability density function (80 bins) of the retrieved abundance of water 
ice, being the value of the model in input equal to 5%. The black line shows 
the resulting distribution in the case is defined according to Eq. 14.11 The 
clear line shows the resulting distribution in the case the a in Eq. 14.11 is not 
coming from the S/N simulator but the observed data is used as weight. The 
theoretical gaussian distributions over plotted (dashed line) have the same mean 
and dispersion of the resulting distributions. The position of the bars are in 
the points of the theoretical distributions such that they contain the 95% of 
the probability distribution. The distance of the bars from the mean value is 
indicated in the figure. 

For real data we could have the situation in which the of the 
S/N calculator underestimates the real noise of the spectrum. This 
is enlighten by an high value (>> 1) of the Xr coming from the 
retrieval procedure. In this case to correctly evaluate the error bars 
of the retrieved parameters, we should scale the calculated error bars 
by the square root of the Xr- 

Other sources of error could add uncertainty on the results, like 
errors on the calibration. Those typically add a hetitious offset to 
the whole spectrum or in some parts of this. Any systematic error 
should be removed before starting the retrieval procedure. Other 
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unpredictable errors should be included in this procedure. As an 
example if we have an estimated uncertainty of 5% on the absolute 
calibration, each synthetic spectrum should include a random offset 
having a a of this amount. 

4.2.2 Errors from geometry uncertainties 

The determination of the angles (incidence and emission) of an area 
subtended by a spatial pixel can presents uncertainty and this should 
not be neglected because they are required by the model (Eq. 13.1|) . 
The shape of the surface is modeled by many facets that have their 
own orientations. The area subtended by a pixel is always larger 
than the area of the facets, and therefore, there are several facets 
for each pixel. 

There are two sources of error that contribute to the uncertainty: 

• error in the estimation of the average angle. This is usually 
performed taking into account the four facets at the corners 
and at the center of the pixel; 

• effective angles could be different with respect to the average 
angles because of the effect of the topographic roughness. This 
effect is predicted by the Hapke model and should not be in¬ 
cluded in this evaluation. 

To predict the effect of the error on geometry we have to calculate 
the probability density plot as described in section 14.2.11 including a 
random variation in the orientation of the surface. 

For each simulation, the normal to the surface is tilted by a ran¬ 
dom angle 6 , that is normally distributed with a hxed standard 
deviation ae. The direction of the tilt is determined by the (fi angle 
that is uniformly distributed. The cosines of the angles required by 
the Hapke reflectance model are calculated following the spherical 
geometry as follows: 

To TILT = To X cos{9) -f sin{i) x sin{6) x cos((/?i) (4.2) 

Rtilt = /i X cos{9) + sin{e) x sin{9) x cos((pe) (4.3) 


where: 


T>i = 

ipi = 27r - (p. 


ipe = (p’ - (f, for (p <(p' 

Pe = p- p', for p> p' 
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ip = arcos {[cos{g) — rq x R\/[sin{i) x sm(e)]} 

being p and p' respectively the azimuth of the tilt and the az¬ 
imuth of the emission direction in the reference frame whose z axis 
is the normal to the surface, that are counted counterclockwise with 
the azimuth of the incidence direction equal to 0° (see Figure 03]). 



Figure 4.4: Schematic diagram of tilt coordinates on an element of surface S. 
Being N the normal to the surface, Nt the normal tilted, i and e the incidence 
and emission angle, ix and ex the incidence and emission angle after the tilt of 
the normal, y and ip' respectively the azimuth of the tilt and the azimuth of the 
emission direction from the direction of the incidence angle; 9 is the tilt angle, 
and g is the phase angle. 

The procedure discussed in section 14.2.11 is performed including 
a random tilt of the surface with a fixed erg, in a way to modify the 
value of the incidence and emission angles, and therefore the cosines 
required by the Hapke reflectance model. 

The dispersion of the retrieved parameters is thus a measure for 
errors that can be committed in the determination of these param¬ 
eters, due to both the signal error and the observation geometry 
uncertainties. 

We can include the evaluation of the error due to geometry infor¬ 
mation if all the spectra are coming from the same surface region. 
In case we consider spectra from different surface regions, in order 
to perform a phase function £t (assuming the same compositional 
properties), the uncertainties on the geometry have different effects 
for each spectrum. In this case this last tool is no longer applicable 
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and the errors on geometry can be considered as a farther sonrce 
of noise that increases the a in Eq. 14.11 We shonld then apply the 
procednre discnssed in section 14.2.11 and consider the valne of the 
Xji coming from the best £t to correctly evalnate the error bars of 
the retrieved parameters. This is the case of onr Lntetia analysis, 
discnssed in chapter 16.11 
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Chapter 5 

Deep Impact-HRI data 
analysis 


On 2005 July 4, the Deep Impact mission carried out an active 
planetary experiment by delivering a high-speed impactor against 
the nucleus of Comet 9P/Tempell. The mission excavated mate¬ 
rial from well below the surface, making it available for study both 
by the Deep Impact flyby spacecraft and by remote sensing from 
Earth, both on the ground and in Earth orbit, and from elsewhere 
in the Solar System (Roset ta, Spitzer Space Observato ry, Solar and 
Heliospheric Observatory) ( A’Hearn and Combi . 20071) . 

The primary goal of Deep Impact was to understand the differ¬ 
ence between the surface of a cometary nucleus and its interior and 
thus to understand how the material released spontaneously by the 
nucleus is related to the primitive volatiles that were present in the 
protoplanetary disk at the age of formation. The other important 
goal was to understand the phys i cal pr operties of the outer layers of 
the comet flA’Hearn and Combi 120071) . 

9P/Tempel 1 is a Jupiter-family comet discovered in 1867 by 
Ernst Wilhelm Leberecht Tempel. With an orbital period at the 
time of 5.68 years, it was observed on two subsequent apparitions. 
However, in 1881 owing to a 0.55 AU pass by Jupiter, its orbital 
period was increased to 6.5 years and its peri helion rai s ed fro m 1.77 
to 2.07 AU. The comet remained “lost” until Marsden ( 1963 ) calcu¬ 
lated that close approaches of the comet to Jupiter in 1941 and 1953 
should have decreased the perihelion distance and orbital period. As 
a result the comet was recovered by Roemer in 1967 and has been 
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observed all of its subsequent perihelion apparitions. Tempel 1 now 
has an orbital period of 5.5 yea rs and a perihelion distance of 1.5 
AU f A’Hearn and Combil . 200?! ). 

Deep Impact was launched from Cape Canaveral Air Force Sta¬ 
tion on 2005 January 12 and it took a short, six-month trajectory 
to position itself to be overtaken by Comet Tempel 1 on the other 
side of the Solar System on 2005 July 4. The encounter occurred one 
day before the comet reached perihelion and 3 days before the comet 
reached its desce nding node, thus opt i miziii g the energy delivered 
in the encounter flA’ Hearn and Combil . 120071) . 

The spacecraft consists of two main sections, the 370-kg copper- 
core “Smart Impactor” that impacted the comet, and the “Flyby” 
section, which imaged the comet from a safe distance during the 
encounter with Tempel 1. 


Flyby 
spacecraft 

Medium 
resolution 
imager 


High-gain antenna 
with 2-axis gimbal 


7.5 m^ 
solar ^ray 



/ 

High 

resolution 

imager 


Smart, 
instrumented impactor 


Figure 5.1: Deep Impact spacecraft. 

The Flyby spacecraft (see Figure 15.ip is about 3.2 meters long, 
1.7 meters wide and 2.3 meters high. It includes two solar pan¬ 
els, a debris shield, and several science instruments for imaging, 
infrared spectroscopy, and optical navigation to its destination near 
the comet. The spacecraft also carried two cameras, the High Res¬ 
olution Imager (HRI), and the Medium Resolution Imager (MRI). 
The HRI is an imaging device that combines a visible-light camera 
with a hlter wheel, and an imaging infrared spectrometer that oper- 
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ates on a spectral range from 1.05 to 4.8 pm. It has been optimized 
for observing the comet’s nucleus. It consisted of a long-focal-length 
telescope followed by a dichroic beam splitter that reflected (0.3 to 
1.0 microns) light through a hlter wheel to a CCD for direct, optical 
imaging (HRIV). The beam splitter transmits the near infr ared light 
to a 2-prisms spectrometer (HRII) flHampton et ahl. 1200511 . 


Table 5.1: HRII Spatial characteristics 


Physical Pixel Size 

36 micrometers 

Effective Pixel FOV 

10.0 microradians 

Effective FOV 

2.5 milliradians or 0.15° 


Table 5.2: HRII Spectral characteristics 


Effective FOV 

10.0 microradians (slitwidth) 

Spectral range 

1.05 to 4.8 fj,m 

Resolution 

R = 200 - 800 


Table 5.3: HRIV 


Pixel Size 

21 micrometers 

Pixel FOV 

2.0 microradians 

Instrument FOV 

2.0 milliradians or 0.118° 

Surface Scale 

1.4 meters/pixel at 700 km 


The MRI is the backup device, and was used primarily for nav¬ 
igation during the hnal 10-day approach. It also has a hlter wheel, 
with a slightly different set of hlters. 

The Impactor section of the spacecraft contains an instrument 
that is optically identical to the MRI, called the Impactor Targeting 
Sensor (ITS), but without the hlter wheel. Its dual purpose was to 
sense the Impactor’s trajectory, which could then be adjusted up 
to four times between release and impact, and to image the comet 
from close range. 

As the Impactor neared the comet’s surface, this camera took 
high-resolution pictures of the nucleus (as good as 0.2 meters per 
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pixel) that were transmitted in real-time to the Flyby spacecraft 
before it and the Impactor were destroyed 

(http://en.wikipedia.org/wiki/Deep_Impact_(spacecraft)). 

The observations of the impact show distinct signs of layering 
within the excavated material and provide an npper limit to the 
strength in the snrface layer showing that the material is mnch 
weaker than ice and that the nnclens is highly porous. 

The appearance of ice in the coma immediately after the impact, 
requires that ice must be present in the hrst outer layers of the 
surface at the impact site, a nd it is important to note that fairly 
pure ice grains were present flA’ Hearn and Combil . 120071) . 

Key results unrelated to the impact include the determination 
of a very low thermal inertia, a very high frequency of spontaneous 
outbursts, apparently correlated with rotational phase in the sev¬ 
eral weeks preceding the impact, and the ex istence of small patches 
of ice on the surface ( Sunshine et ahl. l2006l ) t hat cannot be the pri¬ 
mary source of the water observed in the coma ( A’Hearn and Combil . 

2ffi3). 


Although the impactor spacecraft was destroyed, the flyby space¬ 
craft and its instruments remained healthy in its 3-year, heliocentric 
orbit after completion of the mission. After that the Deep Impact 
flyby spacecraft was retargeted to comet 103P/Hartley 2 as part of 
an extended mission named EPOXI (Extrasolar Planet Observation 
and Deep Impact Extended Investigation). 

The closest approach to Hartley 2 at 694 km happened on 2010 
November 4th, 1 week after perihelion passage and at 1.064 astro¬ 
nomical units (AU) from the Sun. Observations of the comet were 
carried out for 2 months on approach (5 September to 4 November) 
and for 3 weeks on departure (4 to 26 November) , during which 


more than 105 images and spectra were obtained (lA’Hearn et ah 


20111 ). 


Prior remote sensing showed that Hartley 2’s nucleus has an aver¬ 
age radius 1/5 that of comet Tempel I’s nucleus, yet it releases more 
gas per unit time at perihelion, even when allowing for the smaller 
perihelion distance of Hartley 2 (1.059 versus 1.506 AU). This puts 
Hartley 2 in a different class of activity than Tempel 1 or any of 
the other comets visited by spacecraft. The two comets have very 
different surface topography (Fig. 15.21) . but whether this is related 
to the hyperactivity or other processes is still being investigated 
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Figure 5.2: Comparison of a small part of (left) Tempel 1 with (right) Hartley 
2 at approximately the same image scale and with nearly identical instruments. 
Left: Impactor Targeting Sensor (ITS) image iv900067 5, 9.1 m/pixel. Right : 


MRI image mv5004032, 8.5 m/pixel. Sun is to the right (jA’Hearn et al.Ll2011r) . 


(lA’Hearn et al.l. 120111) . 


5.1 Tempel 1 


5.1.1 Spectral analysis 


The calibrated radiance spectra are available through the Plane¬ 
tary Data System website (http://pds.jpl.nasa.gov/). During the 
flyby of the comet many hyperspectral cubes were taken by the 
HRI instrument with increasing resolution. We have analyzed the 
most resolved cube which present the whole nucleus: Exposure ID = 
9000036, with a resolution of 164 - 152 m/px (geometry information 
on this acquisition are courtesy of Tony Farnham). After the impact 
other cubes with a maximum resolution of 16 m/px were acquired 
but they show only a little portion of the surface, contaminated by 
the impact ejecta and without geometry information available. So 
they are discarded for the following spectrophotometric analysis. 

Only data of the spectrometer in the range 1.2 - 4.5 um are 
used, beca use of an higher radiometric accuracy flKlaasen et ah 
2908112ni8[) . 


The first step to perform the analysis is the modeling of the 
radiance in terms of thermal emission and reflected sunlight (see 
Eq. 13.4p . This is useful for both isolating reflectance spectrum 
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and thermal emission. In this section the former contribntion is 
analyzed, in next section the latter. 

The procednre performs the following steps: 

1. It divides the measnred radiance by the solar irradiance cor¬ 
rected for the heliocentric distance {J/D‘^). The result is the 
reflectance added to the thermal emission. 

2. At shorter wavelengths (1.2 - 2.5 pm) the spectrum is not con¬ 
taminated by thermal emission. A linear £t is performed in this 
range to extend the reflectance spectrum at longer wavelength. 
This is reasonable for Tempel 1 spectra that present a quite flat 
trend. Thus for each pixel a spectral reddening is calculated as 
the slope of the reflectance spectrum (see Fig. 15.9p . 

3. The contribution of the reflected sunlight in term of radiance 
is calculated by multiplying the solar irradiance for the inter¬ 
polated reflectance spectrum. 

4. The thermal emission is isolated by the total radiance by sub¬ 
tracting the calculated reflected sunlight contribution (previous 
point) to the measured radiance. 

In this way the two contributions are isolated and are available 
for further analysis (see Fig. 15.3p . 

Isolating the thermal emission is also useful to extend the re¬ 
flectance spectra to longer wavelengths: the thermal emission is 
modeled as discussed in next section and is subtracted to the mea¬ 
sured radiance spectra before dividing it by the solar irradiance. 
In this way the resulting reflectance spectra are free from thermal 
emission signal. 

The data of the camera (HRIV) are available for seven differ¬ 
ent bands from 0.35 to 0.95 pm (see Tab 15.4p . These images are 
rescaled to the same spatial resolution of the hyper-spectral cubes 
obtained by the spectrometer. In this way we are able to extend 
the spectra of each pixel to the visible wavelengths. This task is 
straightforward because all the used data were acquired in a short 
time interval compared with the comet’s rotation period. There¬ 
fore, the dataset analyzed contains only observations taken on the 
same nucleus hemisphere. Since the radiometric accuracy of the 
VIS channel is higher than the IR, the latter has been scaled to the 
former before bridging the spectra of each pixel. This turned to be 
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Figure 5.3: The measured radiance (black line) is modeled in terms of reflected 
sunlight and thermal emission. 



affect the HRII acquisition because of the low spatial resolution of 
the HRII channel, while can be clearly seen for HRIV acquisition 
(see Fig. I5.4p . However, as stated, the HRIV images are rescaled to 
the same spatial resolution of the HRII channel. Thus this effect do 
not affect the analysis. 

The VIS-IR hyperspectral cubes are powerful tools to identify 
icy regions (see Fig. I5.5p . After having assigned thresholds for the 
reflectance at visible wavelengths and for the 2.0 pm band depth is 
possible to identify the ice-rich unit s. In this way we have obtained 
a similar map to that obtained by Sunshine et ah f 2006h (see Fig. 
EH). 


To dehne a SSA of the Dark Terrain we need to exclude the icy 
units from the analysis. Taking into account the Hapke model (Eq. 
13.ip . we can derive a SSA from the average spectrum of reflectance 
of the non-icy units. However, because the available observations 
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Table 5.4: The selected images taken from the HRIV instrument. 


Exposure ID 

Filter (nm) 

Pixel scale (m/pixel) 

Distance (km) 

9000906 

350 

17.95 

8977 

9000908 

450 

17.72 

8861 

9000907 

550 

17.84 

8923 

9000901 

650 

18.67 

9337 

9000903 

750 

18.39 

9196 

9000902 

850 

18.55 

9274 

9000905 

950 

18.14 

9071 



Figure 5.4: Left panel: HRIV acquisition, Exposure ID = 9000901, right panel: 
same image after lowering the spatial resolution to match that of the HRII 
channel. 

are all taken at a single phase angle (63°) we have to assume the 
value of the parameters in the model. At this specihc phase angle 
we assumed that the opposition effect (signihcant at small phase 
angles) and the shadowing caused by the roughness (signihcant at 
high phase angles) are negligible. Moreover, we hx the porosity 
parameter and the single particle phase function to 1. In this way 
we have derived the SSA by inverting the Eq. 13.11 and we have 
performed a linear interpolation of it (see Fig. 13.ip 

The extrapolation of the SSA (ca) is useful in modeling the re- 
hectance spectra as a mixture of Dark Terrain and pure water ice. 
The SSA of pure water ice is obtained from optical constants by 
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Figure 5.5: Average spectra of nucleus regions shown in Fig. 15.61 



Figure 5.6: Panel A: image of the nucleus of the comet at A = 0.35/rTO (exposure 
ID 9000906), rescaled to the same resolution of the hyperspectral cubes. Panel 
B: spectral classification of the comet’s nucleus: in blue the icy regions, in green 
the low reddening regions, in orange the high reddening regions, in grey the 
regions which presents unfavourable angles for the analysis. 


Clark et al.l (120121) , iMastrapa et al.l (120081) and IWarrenI (119841) , cov¬ 


ering a large range of temperatnre (40 - 260 K), and different states 
(crystalline and amorphons). In the final fit there are not remark¬ 
able differences in the resnlts nsing different optical constants. As 
shown in Fig. 15.71 the low amonnt of surface ice (and the instru¬ 
mental sensitivity) don’t allow us to infer the water ice state. In 
particular the 1.65 pm secondary absorption feature, characteristic 
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of the crystalline form is not evident on Tempel 1 spectra. 




Figure 5.7: Intimate (left) and areal (right) mixing results: icy region (corre¬ 
sponding to pixels marked in blue in Fig. 15.61) average observed spectrum (black 
curve); areal mixing spectra with 2.4, 1.2 and 0.4% of water ice are indicated by 
red, green and blue curves, respectively. The best-fit are obtained respectively 
with 1.2% and 0.4% water ice-rich spectrum. 

Adopting the areal mixing we obtain from the spectral hts an 
average grain size of 30 ± 20 pm in diameter with a percentage of 
water ice of 0.4 ± 0.2%. In the case of intimate mixing we obtain a 
grain size of 70±40/im with a percentage of water ice of 1.0 ±0.5%. 
The spectra corresponding to the icy regions are well-htted with 
both mixing methods. 


Water ice abundance (%) 



sample 


Figure 5.8: Water ice abundance map, in case we model the mixing between 
Dark Terrain and pure water ice to be intimate. 
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As a result of the mixing modeling we obtain a map of water ice 
abundance (see Fig. 15.81) . whose abundance slightly differs if we take 
into account an areal or intimate mixing (see Fig. 15.7p . An intimate 
mixing means that the ice and non-icy components are thermally 
coupled. This has implications on ice sublimation as discussed in 
section 15.1.31 


5.1.2 Thermal emission analysis 


The dehnition of the single scattering albedo (u) discussed in pre¬ 
vious section is also useful to constrain the thermal properties of 
the material making up the comet surface. In particular here we 
are interested to deriy e the integrated emissivity {eh)- Following 


Davidsson et ahl (120091 1 it can be approximated with: 


et^l-Ah (5.1) 

where Ai, is the bolometric albedo dehned as the average of the 
spectral spherical albedo As(A) weighted by the spectral irradiance 
of the Sun Js(A): 


Ah = 


r J.(A)<iA 


( 6 . 2 ) 


and the spherical albedo (A^) is the total fraction of incident irradi¬ 
ance scattered by a body i nto all directions. It can be approximated 
by (Eq. 11.45b in Hapke (2012b)): 


As 


1-7 

1-^7 


1 - 


1 7 


31 + 7 


(5.3) 


where 7 = \/l — uj. 

Although the integral to derive Ah should be calculated from 0 
to inhnity, most of the solar radiation comes from the wavelengths 
covered by the channels HRIV and HRII (from ~ 0.2 to ~ 5pm ), 
hence an integration on this range is a good approximation. 

The overall calculation brings to the value of eh = 0.987 

To calculate the Te niperature and th e emis sivity from the data it 
is followed the w ork ofiDavidsson et ahl (120091 ). With respect to the 
previous work of Groussin et ahl ( 2007 ) he admitted the possibility 
of a lower emissivity than that resulting from the very dark surface 
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of the comet, which should be close to 1. The lower emissivity 
is justihed by the effect of roughness which produces self heating 
effects (see Eq. 13.81) which are directly linked to the roughness of 
the surface. 

Following the method discussed in section 14.11 a best £t is per¬ 
formed with the model of thermal emission (Planck’s law) to retrieve 
effective emissivity (ce//) and Temperature (T): 


where: 


B 


^e// ^ 


2hc^ 1 

(A ■ 10 gA-io-e kg T) 


X 10“® 


(5.4) 


the factor 10 ® is needed for a A in input in [p,m] and the 
dimensions in output to be: \W/sr//m^] 

h = 6.62 ■ 10“^^ J • s is the Planck constant 

Kb = 1.38 • 10“^^ J/iP is the Boltzmann constant 

c = 2.998 • 10® m/s is the speed of light 


The effective emissiv ity Cg// is linked to the self heating (,^) by 
f Pavidsson et ah . 20091 ): 




(5.5) 


Thanks to the retrieval for each spatial pixel, maps of tempera¬ 
ture, effective emissivity and self heating are obtained. 

The maps ob t ained (see Fig. 15.9p are very similar to those in 
Davidsson et ah ( 2009 1. 

The position of the ice-rich units (see Fig. 15.61 and Fig. 15.8p is 
coincident with that of the lowest reddening and lowest Temperature 
regions (see Fig. 15.9p . confirming the consistency of the maps. 


5.1.3 Thermal emission analysis: alternative approach 

The interpretation of the results of thermal analysis exposed in the 
pr evious section poses so me difhculties. 


Sunshine et al.l (j2006l ) pointed out that the presence of water ice 


would imply temperature lower than the free sublimation temper¬ 
ature (~ 200K), while resulting temperatures are close to 300 K. 
A possible explanation is that water ice is physically and thermally 
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Figure 5.9: Upper left panel: temperature T{K); upper right panel: effective 
emissivity £eff', lower right panel: self heating; lower left panel: reddening 
(%/kA), it is calculated with respect the value of the reflectance at A = 1.7 ytm 
for a direct comparison with the work of IPavidsson et riL ( 20091 ). 


decoupled from the rest of the surface at higher temperature, and 
the water ice involves only a little portion of the surface subtended 
by the pixels. 


Another difficulty rises when we try to link the self heating pa¬ 
rameter (,^) to the mean slope (9), being both related to the rough¬ 
ness of the s urface. The two parameters can be linke d univocally as 
propo sed by LagerrosI ( 1997 1. However, as stressed by Davidsson et al.l 
f 2009[l . from the resulting mean th e mean 9 shoul d be 55°, signif¬ 
icantly above the value obtained by iLi et ahl (120061 ): 9 = 16°, and 
well beyond the assumptio n of lit t le angle s for the derivation of the 
mean slope parameter bv iHapkel fj2012bl) . Possible explanation of 
the discrepancy between the two parameters is that they measure 
the roughness at fundamentally different size scales: several works 
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indicate that the mean slope parameter (6) measnres ronghness at 
the smallest size scales where the shadows still exist, namely snb- 
millimeter scale, while the self heating parameter (^) conld measnre 
ronghness where this effect is important, that is decimeter scale or 
larger. 

Another possible explanation is that the model neglects the ef¬ 
fect of the shadowing caused by the roughness of the surface. To 
enlighten this point we can make a meaningful example: let’s take 
the case in which the half of the surface with an emissivity e ~ 1 
subtended by a pixel is at temperature < 160 K because it is in 
shadow. Below this temperature, a spectrometer working at wave¬ 
lengths up to 5 pm loses sensitivity to the thermal emission. Let’s 
suppose the other half of the surface subtended by the pixel is at a 
temperature of 300 K. Let’s also suppose the self heating effects are 
negligible. The radiance emitting from such a pixel, in the range 1 - 
5 pm would be equal to the half of the radiance of a black body at 
300 K. The temperature retrieved by the model exposed in previous 
section would be 300 K, with an effective emissivity of 0.5, even if 
self heating effects are negligible. This effect could be enhanced by 
observing at higher phase angle, where shadow effects prevail. In 
this case we could overestimate the mean temperature and the mean 
roughness of the surface region subtended by a pixel. 

Here a different approach in the interpretation of the thermal 
emission is presented. It can overcome most of the problems just 
discussed. The intent of this approach is to admit the surface sub¬ 
tended by a pixel to be made by different portions emitting at dif¬ 
ferent temperatures because of the roughness effects. The thermal 
emission is modeled as a sum of Planck functions at different tem¬ 
peratures in the range T = 0 - 400 K (with an integrated emissivity 
Ih as calculated in section l5.1.2p . weighted by a Gaussian function 
G{T), whose mean and standard deviation are the free parameters 
to retrieve. The mean value represents the most frequent value of 
temperature on the surface subtended by the pixel, while the stan¬ 
dard deviation is the measure of the dispersion of the temperature 
values caused by roughness effects. Results are mapped in hgure 

Eini 

The map in hgure 15.101 (left panel) shows temperatures signih- 
cantly lower that those in hgure ED (upper left panel). However, it 
should be stressed that the results coming from the two models are 
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Figure 5.10: Mean (left) and standard deviation (right) map of the temperature 
on the surface, according to the Gaussian distribution of Planck functions model. 


not contradictory because are based on different assumptions. If we 
use the model presented in previous section we should keep in mind 
that the radiance is more affected by thermal emission coming from 
higher temperature regions. This last method is recommended if the 
risk of contamination by shadows is high, and we are interested in 
modeling a distribution of temperatures on scale smaller than that 
of the pixel. 

The icy regions mapped in hgureES present a mean temperature 
T ~ 250, with a standard deviation of dr ~ 30. According to a 
Gaussian distribution the probability to have temperature < 200iG 
with this parameters is ~ 5%. Thus, such a portion of the surface 
is compatible with the presence of the exposed water ice without 
risk of sublimation. This result is consistent with the abundance of 
water ice retrieved in those regions (< 2%). 


5.2 Hartley 2 

5.2.1 Spectral analysis 

The calibrated radiance spectra are available through the Planetary 
Data System website (http://pds.jpl.nasa.gov/). During the flyby 
of the comet many hyperspectral cubes were taken by the HRI in¬ 
strument. As in the case of Tempel 1 analysis, we take into account 
the most resolved cube which present the whole nucleus: Exposure 
ID = 5005001, with a resolution of 28.5 m/px (geometry information 
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on this acquisition are courtesy of Tony Farnham). 

The radiance is modeled in terms of thermal emission and re¬ 
flected sunlight (see Eq. I3.4p following the steps discussed in sec¬ 
tion ISHITl The spectra of Hartley 2 appear to be more differentiated 
across the surface than those of Tempel 1. For this reason a different 
approach for the spectral analysis is followed: the spe ctra are mod¬ 


eled a s binary intimate mi xtures of amorphous carbon fIZnbko et ah 


1996) and pu re water ice f Clark et ah . 20121: Mastrapa et ah . 2008 


Warrenl. 1984 ). thus avoiding the dehnition of a dark terrain, as done 
for Tempel 1. 

Following the method discussed in section 14.11 the abundance 
and grain size of the two endmembers are retrieved for each pixel of 
the comet’s nucleus. 


Water ice abundance (%) Water ice grain size (gm) 



Figure 5.11: Water ice abundance and grain size map of Hartley 2, in case we 
model the mixing between amorphous carbon and pure water ice to be intimate. 


The abundances retrieved can be intended as the area covered by 
that material as usual, but because the grain size of the two end- 
members are independent each other, it cannot be also intended as 
a numerical density anymore. 

Figure 15.111 shows large regions on the nucleus surface where 
the water ice is relatively abundant (up to ~ 20%). The position 
of these regions on the nucleus surface match with the distribu¬ 
tion of water ice in the inne r coma as a product of activity of the 
comet flA’ Hearn et al.l . 120111 1. Moreover, the region where exposed 
ice present larger grain size seems to be correlated with the most 
abundant and hetero geneous {H 2 O, CO 2 , CN) ambient outgassing 


f A’Hearn et ah . 20 111) 
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Figure 5.12: Upper panel: the measured spectrum and the retrieved model. 
Lower panel: residuals analysis. The spectrum of residual is smoothed by a 
running box. A signifi cant residual f> 2a) at A ^ 2.7 /rm is related with the 
water vapour emission ( A’Hearn et al. . 201lh which is not taken into account in 
the model. The best fit is performed skipping the spectral range 2.55 - 2.8 /rm 
(dashed spectrum), related to this water vapour emission. Retrieved water ice 
abundance: 1 %; retrieved grain sizes: water ice = 100 /im, amorphous carbon 
= 4 fim. 


Figure |5ll2] (upper panel) shows the measured spectrum and the 
best £t model, from which abundance and grain size of water ice 
are retrieved. The spectrum is related to the clearest region of right 
panel of hgure 15.111 which present larger grain size of the water ice 
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Figure 5.13: Upper panel: the measured spectrum and the retrieved model. 
Lower panel: residuals analysis. The spectrum of residual is smoothed by a 
running box. A signifi cant residual (> 2a) at A ^ 2.7 fim is related with the 
water vapour emission ( A’Hearn et al. , 20 111) which is not taken into account in 
the model. The best fit is performed skipping the spectral range 2.55 - 2.8 /rm 
(dashed spectrum), related to this water vapour emission. Retrieved water ice 
abundance: 8 %; retrieved grain sizes: water ice = 7 /xm, amorphous carbon = 
2 fim. 


particles. The analysis of the goodness of the model is given in the 
lower panel. 

Figure 15.131 shows the same analysis of fignre [5T2] for a pixel in 
a dark region of right panel of fignre 15.111 (small grain size region). 
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The difference of the spectra coming from the two regions (Fig. 15.121 
and 15.13]) is clearly visible: the former presents the three absorption 
bands of water ice (A ~ 1.5, 2.0 and 3.0 /rm), the latter present a 
large absorption band at A ~ 3.0 fim, being the other two weak or 
absent. 


5.2.2 Thermal emission analysis 

As discussed in sec tion 15.1.21 the th e rmal analysis is performed fol¬ 
lowing the work of Davidsson et ah ( 20091) . The spectrum of single 
scattering albedo for comet Harley 2 is extracted from reflectance 
spectra in the same way of Tempel 1 (see section l5.1.1|) . The slope 
and the level of the albedo (cn) of the two comets are very similar 
as shown in Fig. 15.141 For this reason we assume for Hartley 2 
the same integrated emissivity Ih already calculated for Tempel 1 in 
section 15.1.21 



Figure 5.14: The single scattering albedo (w) calculated for the two comets. The 
high level of uncertainty coming from the thermal emission removal prevents the 
calculation of the SSA beyond 3.5 /rm. The SSA derived from Hartley 2 data 
presents contamination of water vapour emission (A ~ 2.7) not removed here. 
Superimposed is the linear interpolation of the Tempel 1 SSA, which has been 
used for the spectra simulation in Chapter [3| The two slopes are related to the 
two different channels of the instrument: HRIV and HRII. 
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Maps of Temperature, effective emissivity, self heating and red¬ 
dening are calculated as discussed in section 15.1.21 and are shown in 
Fig. [5T51 


Temperature (K) 


Effective Emissivity 



sample 


Reddening (%/ kA'^) 


Seif Heating 



sample 



sample 



Figure 5.15: Upper left panel: temperature T[K)] upper right panel: effective 
emissivity Ce//; lower right panel: self heating; lower left panel: reddening 
(%/kA), it is calculated with respect the value of the reflectance at A = 1.7 
for a direct comparison with the work of lPavidsson et ^ ( 20091) . 


The reddening map (lower left panel) shows results in the same 
range than those of Tempel 1 (see Fig. 15.151) . The two horizontal 
stripes are products of artifacts on the spectra. 

As in the case of Tempel 1 the map of temperature (upper left 
panel) is consistent with the direction of the Sun. The map of ef¬ 
fective emissivity (upper right panel) show very low values, that are 
even more problematic in their interpretation than those of Tempel 
1. The related map of self heating shows as a consequence very high 
values. Moreover, higher values are located in regions weakly illu¬ 
minated by the Sun, showing that shadowing effects are important 
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and produce misleading results. 

The alternative approach discussed in section 15.1.31 is applied to 
Hartley 2 data. Results are shown in Fig. 15.161 



Figure 5.16: Mean (left) and standard deviation (right) map of the temperature 
on the surface, according to the Gaussian distribution of Planck functions model 
(see section H.1.3l l 

As in the case of Tempel 1, temperature map in Fig. 15.161 (left 
panel) shows lower values than those in Fig. 15.151 fupper left panel), 
because of the different modeling. However, the alternative ap¬ 
proach can better explain the presence of exposed water ice up to 
20 % in abundance (see Fig. 15.lip , having those regions a mean 
temperature < 200 K. 

The dispersion of Temperature distribution (right panel) shows 
a mean value ctt ~ 60. This higher value with respect to that re¬ 
sulting from Tempel 1 {ax ~ 30) could be mostly related to the 
higher phase angle during Hartley 2 data acquisition (84" against 
63°), which produce more shadows on the surface and so more tem¬ 
perature differentiation. 
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Chapter 6 

Rosetta-VIRTIS data 
analysis 


6.1 Lutetia photometric analysis 


On 10th of July 2010 the Rosetta spacecraft was directed to a fly¬ 
by with the main belt asteroid 21 Lutetia, a large sized asteroid of 
~ 100 km in diameter. Its ancient surface age (determined from 
crater counting) coupled with its complex geology and high den- 
sity suggests that L utetia is most likely a primordial planetesimal 


fjS jerks et al.l . l201lh . This was also conhrmed by t he spectroscopic 


obser vations performed by the VIRTIS instrument flOoradini et ah 


201 ll ). which have shown no absorption features, of either silicates 
or hydrated minerals (see Figure lOTl) . 


The availability of accurate geometric information (iKeihm et ah 


2 OI 2 I ) (incidence, emission and phase angles) for each pixel in the 
VIRTIS acquisitions and in all phases of the fly-by gives us the 
opportunity to perform a disk-resolved photometric analysis. 


Analysis of the normalized spectral variation across the observed 
surface at highest resolution shows a remarkable uniformity of the 
surface spectral properties. For this reason we can assume all the 
surface having the same photometric characteristics, and we can use 
the large amount of available spectra in different viewing geometry 
to improve the retrieval of the photometric parameters. 
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Figure 6.1: Lutetia mean reflectance spectrum in the 1.0 - 5.0 fim spectral range. 
Where the line is thin there are instrumental artifacts in the spectrum; the rise 
above 4 fim is due to the surface thermal emission. In light color are shown the 
Icr spectra. The panel inside the plot shows a detail of the reflectance within 
the wavelength range taken into account for the analysis. 


6.1.1 Data and reduction 


During the close encounter with Lutetia, VIRTIS obtained different 
high resolution disk-resolved hyperspectral cubes from IR and VIS 
channels, at a wide range of phase angles. To study the photometric 
properties, we focused on the IR channel from 1.1 to 2.4 pm which 
is a reliable spectral region for the absolute photometric calibration 
(see Figure 16.ip . and do not present contamination from thermal 
emission. The image s were calibrated t hrough the standard VIRTIS 
calibration pipeline ( Filacchionel. 2006 ) to reflectance units. Out of 
statistics pixels are removed by applying the despiking procedure 
discussed in section 16.2.11 

Table 16.11 summarizes the hyperspectral cubes with high spatial 
resolution taken into account for the analysis (see also Fig. 


This dataset spans across a wide range of phase angles 0.01°—136.4°, 
however, as the angular radius of the Sun at the Lutetia distance 
(2.72 AU) is 0.1° observations at smaller phases do not improve the 
knowledge of the photometric properties of the object. Pixels with 
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phase angles smaller than 0.1" are then discarded. 

Table 6.1: The four cubes taken into account for the analysis. They are identified 
with their SpaceCraft Event Time (SCET). Each line of the cubes corresponds 
to a specific phase angle. In the last column the number of spatial pixels selected 
for the analysis. 


SCET 

km/px 

Phase angle range 

Lines 

Selected pixels 

00237396952 

1.7 

12.2°- 136.4° 

175 

9177 

00237396112 

3.5 

0.0°-2.0° 

80 

1545 

00237395857 

5.6 

3.0°-3.2° 

23 

155 

00237394252 

6.2 

4.0° - 7.5° 

296 

620 


The information on the geometry is used for a hrst selection of 
the spatial pixels to be taken into account. Pixels near the limb 
(within 2 pixels) are excluded to avoid pixels not completely hlled 
by the projected area. Moreover, pixels whose incidence or emission 
angle overcome 70" are excluded to avoid the non-linear effect of 
the topographic roughness which causes large uncertainties for the 
model. 

Pixels are further analyzed line by line: for each pixel the average 
of the reflectance over the considered spectral range is taken into ac¬ 
count. This quantity is studied as a function of the geometry accord¬ 
ing to the Pommel-Seeliger law. A linear correlation is expected, be¬ 
cause the phase angle is hxed within the lines: < r >~ [/ro/(/i + /xo)] 
(see Eq. 13.ip . The scatter of the pixels from a linear trend is due 
mostly to geometry information errors. The pixels whose average 
signal overcome 2ct of the scatter are excluded from the analysis. 

6.1.2 Method and models 

Given the flatness of the spectra in the considered range (1.1 — 
2.4 fim) and the need to mitigate the signal uncertainties due to 
the noise, best hts of the phase function are performed taking the 
average over each 10 bands, that is each 0.1 fim. 

The Hapke model (Eq. 13. ip is used for the hrst three hyperspec- 
tral cubes in Table lOTl The least resolved of the four (00237394252) 
presents 14 scans of the full nucleus, permitting a full-disk analysis. 
Although the geometries of this cube are available it is preferred to 
treat it with an integrated formula that takes into account the whole 
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Figure 6.2: The four hyperspectral cubes analyzed at A = 2.0 ^m. They are 
(from up to down) the four listed in table 16.11 


disk, losing the information on the shape model. This permits to 
avoid large uncertainties of geometry information, and to test the 
consistency of the Hapke formula in the full disk version: 
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< FDR > = 


fApdA 

( K/8 |l — sm(^) X tan(^) x /n[cot(^)]| x 
{[1 + BsoBs{g)p{g) - l]w + 4ro(l - tq)} + 

X [sin{g) + (vr - g)cos{g)] ) x (1 + BcoBc{g)) x S{g, 9) / n 

( 6 . 1 ) 


Where: 


ro = (1 - 7 )/(l + 7 ) 

7 = 7/(1 -^) 

S{g,) is the correction for the shadowing function dep ending 
on the roughness parameter 6 (Eq 12.61 pag. 330 in 
fl 2012 b[) h 


Hapkel 


Eg■ 16 ■ 1 1 is straigh tforwardly derived by Eq. 11.35 - 11.42, p. 
299 in iHapkel f 2012bl) and assumes a spherical surface. The integra¬ 
tion gives as a result the reflectance of the full disk averaged over 
the projected surface. This quantity is directly comparable to the 
measured reflectance averaged over all the pixels of the disk. 

The phase function £t is performed taking into account all the 
10877 spectra coming from the first three cubes of Table [ 6 Tl apply¬ 
ing Eq. 13.11 and the 14 spectra resulting from the average full disk 
reflectance coming from the last cube of Table 16.11 using Eq. 16.11 
(see Figure ES])- 

The 9 parameters of the Hapke model hardly fit a unique solu- 
tio n. This problem is recurrent in this kind of analysis as enlighten 
by Helfenstein and Shepardl ( 2011 1. In particular three couple of 
parameters cannot be constrained separately: the porosity parame¬ 
ter K and the single scattering albedo a;, and the four parameters 
concerning the opposition effect {Bsq, hs and Bqn, he). _ 


To constrain K it is followed the suggestion of I Helfenstein and Shepard 
f 201l[l . linking the porosity parameter K, to the angular amplitude 
of the SHOE hs, which also depends on the porosity. It is used 
the Eq. 9.26 pag. 234 in IHapkel fl2012bll . that is derived for equant 
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Figure 6.3: Left panel: reflectance centered in 1.142 ^m, filter 0.1 fim wide, 
plotted against phase angle for all selected spectra. Scatter of the points is due 
to different observation geometry (incidence and emission angle). In the zoom 
is shown a detail for small phase angles, three of the four hyperspectral cubes 
cover this range. From 4° to 7.5° the points are the result of the average over 
the whole disk (14 scans). Right panel: the best fit of the model. It is able 
to account for different observation geometry. Most of the discrepancy between 
data and model are probably due to errors in the geometry information. 


particles, larger than the wavelength, and with a narrow size distri¬ 
bution. 


To constrain the opposition effect parameters it is followed ICiarniello et ah 


f 201ll) by eliminating Bco and he from the fit but still taking into 
account a possible contribution of the CBOE by permitting Bso to 
overcome 1. It is preferred to model SHOE instead of CBOE because 
the typical pre dicted angular width of the SHOE is larger than that 
of the CBOE f Hapkd . 2012bl ). Moreover, the relatively dark surface 
make the opposition effect more affected by SHOE as it is linked to 
single scattering properties, while the CBOE is linked to multiple 
scattering. 


6.1.3 Results and discussion 

The best fit is performed with the method explained in section 14.11 
Figure 16^ shows the best fit for one wavelength. The parameters 
retrieved as a function of the wavelength present a flat trend, as we 
should expect considering the flatness of the spectrum (see Table 
Scatter of the points from this trend are well explained by the 
error bars, that are calculated as explained in section (see Figure 
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EM- 


Table 6.2: Retrieved parameters of the Hapke model as a function of the wave¬ 
length for photometric analysis of Lutetia. In the first line the weighted average 
over the wavelengths range is considered. 


\{pm) 

UJ 

K 

b 

V 

en 

Bo 

hs 

Average 

0.2544 

1.181 

0.263 

0.852 

23.64 

1.575 

0.0525 

1.142 

0.2540 

1.174 

0.276 

0.847 

24.28 

1.598 

0.0490 

1.236 

0.2501 

1.184 

0.265 

0.904 

23.71 

1.667 

0.0539 

1.331 

0.2523 

1.186 

0.264 

0.898 

23.71 

1.640 

0.0547 

1.425 

0.2523 

1.186 

0.259 

0.946 

23.23 

1.629 

0.0544 

1.519 

0.2558 

1.184 

0.265 

0.873 

23.31 

1.602 

0.0536 

1.614 

0.2542 

1.185 

0.260 

0.912 

23.20 

1.604 

0.0544 

1.708 

0.2563 

1.182 

0.263 

0.867 

23.51 

1.582 

0.0527 

1.803 

0.2551 

1.183 

0.263 

0.834 

23.69 

1.580 

0.0531 

1.897 

0.2571 

1.183 

0.260 

0.851 

23.53 

1.590 

0.0533 

1.991 

0.2546 

1.185 

0.256 

0.824 

23.88 

1.634 

0.0544 

2.086 

0.2593 

1.176 

0.265 

0.810 

23.69 

1.508 

0.0499 

2.180 

0.2568 

1.177 

0.264 

0.817 

23.49 

1.490 

0.0505 

2.274 

0.2509 

1.178 

0.260 

0.827 

23.30 

1.492 

0.0509 

2.369 

0.2531 

1.181 

0.259 

0.732 

24.45 

1.515 

0.0524 


Results fT able 16.211 are similar to those comi ng from othe r aster 


oids: St eins fISpiuth et ah . 2012), Annefrank (Hillier et ah . 2011 1. 


Itokawa (iLederer et ah . 2005 : Kitazato et ah . 2008 ). E ros (Clark et ah 


2002 


Domingue et al.l. 120021: iLi et al.l. 120041) . Toutatis (IHudson and Ostrd. 


19981), C astalia (IHudson et al.l. 119971) . Ida (iHelfenstein et al.l. Il 


Gaspra (iHelfenstein et al.l. Il994l) . However, comparisons must be 


made with care because of different channels used for observation, 
and because variations in the models. 

The hlling factor resulting from th e mean retrieve d K parameter 
is: cj) = 0.125 (Eq. 7.45b, pag. 167 in iHapkel f 2012bl) b 

The two parameters of the Heyney-Greenstein function (6, v) 
model the single particle phase function to be mostly backscattering. 
The retrieved value for v (0.852 in average) is always within the 
limits given in Eq. 13.31 (1.50 in average). 

The retrieved Bq overcomes 1, indicating that both SHOE and 
GBOE are present. 

K, 9 and hs should not depend on wavelength, at least in the lit¬ 
tle range considered, because they are linked to the structure of the 
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surface. The next step of the analysis would be to run the retrieval 
by hxing these parameters to their average value, and obtaining the 
other parameters as a function of the wavelength. In this case all the 
values obtained are constant within the errors, because the flatness 
of the spectra, thus this further step does not add any signihcant 
information. 



Wavelength (pm) 


Figure 6.4: The retrieved parameters are plotted as a function of the wave¬ 
length. The mean slope 6 is expressed in radians, like the hs parameter. The b 
parameter is plotted with an offset of -0.1 for the sake of clarity. Error bars are 
calculated as discussed in section IT^ 
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6.2 Comet 67P/Churyumov-Gerasimenko 

6.2.1 Despiking and artifacts removal 

The acquired spectra are always affected by spikes caused by cosmic 
rays and other solar particles hitting the detectors. Their removal 
is specially needed to correctly fit the measured data with models. 
A routine has been developed to address this task. It is devoted to 
deal with spectra showing features at scales much larger than the 
width of the bands of the instrument (~ 10 nm). The procedure 
works as follows: 


1. The users have to set the number of parts in which the spec¬ 
trum has to be divided. Each part will be interpolated with a 
polynomial fit of an order settable by the user. The number 
of parts and the degree of the polynomial fit depends on the 
complexity of the spectrum in terms of spectral features. 

2. The standard deviation of the spectrum with respect to the 
polynomial fit is calculated for each part. The signals overcom¬ 
ing an adjustable number of standard deviations is set equal to 
the polynomial fit. Normally for the spikes removal a number 
of 4 (T is enough. 

3. The previous points are repeated until there are not anymore 
spikes to cut. 


Some parts of the spectra are further aff ected by some art ifacts 
like the crack of the filter in the IR channel (iFilacchiond . l2006h and 
a not proper characterization of the Instrument Transfer Function. 
Many artifacts are proportional to the absolute level of the signal. 
The despiking procedure turned to be useful for artifact removal: 
thanks to the settable level of intensity of despiking, the artifact 
can also be removed. 

However, an intense despiking procedure cannot be applied in 
blind way to all the spectra: the risk to cut real features is high, 
especially if the spectral feature involves few instrumental bands. 

For this reason a map of artifacts is developed for both VIS and 
IR channels. To produce this map the despiking procedure is applied 
to spectral cubes which are assumed to present not spectral features 
at small scale. The map shows the relative difference of the signal 
before and after despiking: 
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Artifact Map {band, sample) 


signal — signal despiked 
signal despiked 


( 6 . 2 ) 



Bands 


Figure 6.5: The map of artifact in the IR channel. It is the median of the Eq. 
16.21 for a large numbers of acquisitions. The crack of the filter is clearly visible. 



Bands 

Figure 6.6: The map of artifact in the VIS channel. It is the median of the Eq. 
16.21 for a large numbers of acquisitions. 

The calculation is performed for each acquisition frame, i.e. each 
line of the cubes, and the median is performed. As a result the maps 
show the artifacts that are signal dependent, thus they can be used 
for their removal: 
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Signal {band, sample) = signal/{I + Artifact Map) (6.3) 

Figure [621 shows the result of the cleaning procedure for a spec¬ 
trum acquired during Lutetia flyby. 

The fact that the maps calculated with Comet 67P CG data 
work well if applied to Lutetia data makes us sure we are removing 
artifacts and not spectral features. 



Figure 6.7: An example of cleaned spectrum of Lutetia after despiking and 
artifact removal. 


6.2.2 Spectral analysis 

A complete spectrophotometric analysis of the nucleus of comet 
67P/C-G is outside the scope of the present work. 

The aim of this section is to show the very hrst results and, as a 
consequence, the hrst problems we have to face. 

The spectra of the resolved nucleus, collected in the hrst phase 
of the mission, are cleaned as discussed in previous section. Sub¬ 
sequently the thermal emission is removed as explained in section 
15.1.11 This allows to isolate the contribution of the rehected solar 
light, namely the rehectance spectra. 

The bridging of the two channels of VIRTIS-M (VIS and IR) 
allows to have a unique view of the spectra from 0.2 to 5 pm, being 
the two spectral channels slightly overlapping at A = 1.0 pm. 
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Despite of the high spatial resolution, the comet’s spectra appear 
to be very uniform along all the surface and basically featureless, 
except for a broad band at 2.9 - 3.7 pm and a tentative band at 0.9 
fim. 

Others distinctive characteristics of the spectra are the very low 
reflectance and two different slopes separated by a knee at about 
1.2 fim. 

Increases of the band depth and width of this broad absorption, 
correlated also to changes in the spectral slopes, occur in active 
areas (neck regions in Fig. 17.ip . 

The method for spectral analysis discussed in section 0] and al¬ 
ready applied for Hartley 2 (see section 15.2. ip has been applied for 
some regions of the comet 67P/C-G, and a typical result is showed 
in Figure ESI 

The spectra are modeled as intimate mi xtures taking i nto ac count 
the optical constants of amorphous carbon fIZubko et al.l.ll996), pur e 


water ice fIClark et al.l. 120121: iMastraoa et al.l. 120081: IWarrenl. 119841) , 


and Titan tholis (llmanaka et al.l . 120121 ). 


The best £t procedure apparently does not need the presence 
of water ice. However, in active areas, where the depth of the ab¬ 
sorption band at 2.9 - 3.7 pm slightly increase, presence of small 
amounts of water ice (1 - 2 %) seems to improve the best £t. 

The retrieved abundance of tholins is ~ 8 %. Even if the slopes 
overall the spectrum are well htted, many features are poorly htted 
as highlighted by the residuals analysis: the model is not able to £t 
the knee at 1.2 fim. In general the measured spectrum is more flat 
than the model in the range 1.2 - 2.5 pm. 

The main issue is represented by the broad band at 2.9 - 3.7 pm: 
it seems the optical constants available to us are not able to correctly 
produce the right position of this band, and its smooth prohle. 

Moreover, from A ~ 4 fim, the signal diverges from the model. 
This could be due again to the not correct optical constants we are 
using to produce the model or, in this case to a not complete removal 
of thermal emission that at these wavelengths provides the same 
contribution of the reflected sunlight to the radiance, and longward 
becomes predominant. 

Finally, at the very little grain sizes retrieved, the approxima¬ 
tion of geometrical optics is not fully satished, and consequently the 
Hapke model could be no longer valid. 
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Figure 6.8: Upper panel: the measured spectrum and the retrieved model. 
Lower panel: residuals analysis. The spectrum of residual is smoothed by a 
running box. Significant residuals (> 2a) affect many parts of the spectrum, 
which is in fact poorly fitted by the model. For a complete discussion refer to 
the text. Retrieved abundances: tholin = 7.5 %, amorphous carbon = 92.5 %; 
retrieved grain sizes: tholin = 2.5 /xm, amorphous carbon = 1.1 fxm. 
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Chapter 7 

Conclusions 


This work has started with a careful analysis of the performance of 
the VIRTIS instrument, based on data acquired during the overall 
cruise phase of Rosetta spacecraft. This analysis is hnalized to a 
better comprehension of the behaviour of the instrument with re¬ 
spect to the signal coming from the comet 67P/CG - main target of 
the mission - under all the observing conditions planned during the 
journey of the comet toward the Sun. 

This has brought to the modeling of the Signal to Noise ratio 
(Chapter 2) which has useful applications for the following analysis. 

To understand how the VIRTIS instrument handles the spectra 
of the comet, “possible” spectra of a comet have been simulated. 
The base for their simulation is the Hapke radiative transfer model 
described in Chapter 3. The modeling of the signal to noise ratio 
has further allowed to add simulated noise to the spectra. 

The tool to produce syntethic noisy spectra is useful expecially in 
the planning phase, to understand under which condition a signal 
is correctly detected by the instrument and therefore to optimize 
the instrument operational parameters depending on the observing 
conditions. Moreover, the analysis of the error on the diagnostic 
absorption bands allow us useful comparisons with real spectra. If a 
certain end-member is not observed we have the opportunity to £x 
an upper limit to its concentration thanks to the method proposed 
in Chapter 3: for a given viewing geometry and instrumental pa¬ 
rameters, the upper limit of its abundance corresponds to the 100% 
error on the band area. 

The modeling of the signal to noise ratio of the instruments is 
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the input of a new method developed to perform spectrophotometric 
analysis. This method, described in Chapter 4 allows an accurate 
and fast retrieval of the model-parameters we are interested in, and 
an estimate of their errors. 

This model has been applied to data obtained from Deep Im¬ 
pact mission during the encounter of the comet Tempel 1 and its 
extended investigation during a second fly-by with the comet Hart¬ 
ley 2. The DI spacecraft carries a spectrometer similar to VIRTIS. 
Results obtained are similar to those obtained by previous works on 
the same targets, showing the reliability of the method. The inves¬ 
tigation of Tempel 1 and Hartley 2 has been fruitful for at least two 
other reasons: new results on the abundances and grains size of wa¬ 
ter ice detected are obtained, and the development of an alternative 
method to perform thermal emission analysis. The latter in partic¬ 
ular allows a better explanation of the presence of exposed water ice 
on the surface of the comets thanks to a modeling of distribution of 
temperatures in the surface subtended by a pixel. 

We have then applied the same methods to the VIRTIS data. 

The photometric model of the Lutetia asteroid have been derived 
applying the developed tools. Results of the photometric parame¬ 
ters retrieved are very similar to recent works on other asteroids, 
using the Hapke’s model. However, the present study differs from 
the previous by the introduction of the K porosity parameter. We 
must emphasi ze that the reliability of this parameter is still under 
investigation. ICiarniello et al.l (1201^ have shown that the updated 
Hapke’s model including K fjHapkd . l2012bl) improves the output at 
medium-large phase angles but also pointed out an inconsistency 
between analytical solution and numerical simulations when oppo¬ 
sition effect and K parameter are modeled at the same time, and the 
latter is >> 1. However, this is not our case because the retrieval 
of porosity K parameter gives back a value close to 1. The porosity 
modelling deserves further investigations, especially because it can 
constrain thermal inertia of the asteroid surface. 

The hnal step of this work coincides with the hrst step of the 
analysis of comet 67P/CG, the main target of the Rosetta mission. 
The hrst result obtained by applying the developed methods to the 
spectra of the comet are very encouraging although still preliminary. 

The very low rehectance of the nucleus as observed by VIRTIS 
coupled with distinct spectral slopes in VIS and IR ranges sug- 
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gests the presence of macromolecular carbon-bearing compounds. 
The very broad absorption in the 2.9 - 3.7 pm is compatible with 
nonvolatile organic macromolecular materials, indicating a complex 
mixture of various types of C-H and/or 0-H chemical groups. 

From the observations performed at heliocentric distance in the 
range 3.6 - 3.3 AU there is no evidence even at the highest spatial res¬ 
olution of 15-25 m/pixel, of ice-rich patches, indicating a generally 
dehydrated nature for the entire surface layer currently illuminated 
by the Sun. 

Further investigations are needed to asses the physical and com¬ 
positional properties of the surface. Attempts should be made to 
understand the origin of its shape in relation to the spectrophoto- 
metric analysis of the three parts in which it can be divided: the 
head, the neck and the body (see Figure ITT]) . 

The time is the fourth dimension in which our analysis will move, 
as well as the comet which will change its heliocentric distance. Al¬ 
though this further complicates the investigation, it could represent 
the key to unveil the secrets of the comet. 
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Figure 7.1: Rosetta navigation camera image taken on 19 

August 2014 at about 79 km from comet 67P/Churyumov- 
Gerasimenko. The comet nucleus is about 4 km across. 

{http://www.esa.int/spaceinimages/Images/2014/08/Comet-on-19-August-2014- 
-NavCam) 
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